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THE COMPRESSIBLE VISCOUS SURFACE-INTERNAL WAVE PROBLEM: 
NONLINEAR RAYLEIGH-TAYLOR INSTABILITY 

JUHI JANG, IAN TICE, AND YANJIN WANG 


Abstract. This paper concerns the dynamics of two layers of compressible, barotropic, viscous 
fluid lying atop one another. The lower fluid is bounded below by a rigid bottom, and the upper 
fluid is bounded above by a trivial fluid of constant pressure. This is a free boundary problem: 
the interfaces between the fluids and above the upper fluid are free to move. The fluids are acted 
on by gravity in the bulk, and at the free interfaces we consider both the case of surface tension 
and the case of no surface forces. We are concerned with the Rayleigh-Taylor instability when 
the upper fluid is heavier than the lower fluid along the equilibrium interface. When the surface 
tension at the free internal interface is below the critical value, we prove that the problem is 
nonlinear unstable. 


1. Introduction 

The Rayleigh-Taylor instability, one of the classic examples of hydrodynamic instability, is 
an interfacial instability between two fluids of different densities that occurs when a heavy fluid 
initially lies above a lighter one in a gravitational field. The instability is well-known since the 
classical work of Rayleigh m and of Taylor BE, and it is one of the fundamental problems 
in fluid dynamics. A general discussion of the physics related to this topic can be found, for 
example, in be- 

The Rayleigh-Taylor instability problem has received a lot of attention in the mathematics 
community due both to its physical importance and to the mathematical challenges it offers. The 
linear Rayleigh-Taylor instability is well understood (see, for instance, Chandrasekhar’s book 
El)- However, there is no general theory that guarantees the passage from linear instability 
to nonlinear instability for PDEs, so the question of nonlinear instability is not immediately 
resolved by the linear analysis. 

In this paper, we are concerned with the nonlinear dynamical Rayleigh-Taylor instability 
of viscous compressible two-fluids having different densities along a free interface, when the 
upper fluid is heavier than the lower fluid along the equilibrium interface. This is the final 
paper in a trio mm that establishes a sharp stability criterion for the compressible viscous 
surface-internal wave problem. 

1.1. Governing equations in Eulerian coordinates. We consider two distinct, immiscible, 
viscous, compressible, barotropic fluids evolving in a moving domain Q(t) = bl + (f) U fl_(f) for 
time t > 0. One fluid (+), called the “upper fluid”, fills the upper domain 

M+(t) = {y € T 2 x R | ri-(y 1 ,y 2 ,t) < y 3 < l + t?+(yi, y 2 , £)}, (1-1) 

and the other fluid (—), called the “lower fluid”, fills the lower domain 

M-(t) = {y G T 2 x R | -b < 2/3 < ri-(yi,y 2 ,t)}. (1.2) 

Here we assume the domains are horizontally periodic by setting T 2 = (27 tLiT) x ( 2nL 2 T ) for 
T = R/Z the usual 1-torus and L\,L 2 > 0 the periodicity lengths. We assume that £,b > 0 
are two fixed and given constants, but the two surface functions rj± are free and unknown. The 
surface r + (t) = {y 3 = £ + rj + (yi, y 2 , £)} is the moving upper boundary of 12+(t) where the upper 
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fluid is in contact with the atmosphere, T_(t) = { 2/3 = ? 7 _ ( 2 / 1 , 2 / 2 , i)} is the moving internal 
interface between the two fluids, and X;, = { 2/3 = — 6 } is the fixed lower boundary of 

The two fluids are described by their density and velocity functions, which are given for each 
t > 0 by p±(-,t) : fi±(t) —> M + and u±(-,t) : f2±(i) —> M 3 , respectively. In each fluid the 
pressure is a function of density: P± = P±(p±) > 0, and the pressure function is assumed to be 
smooth, positive, and strictly increasing. For a vector function u € M 3 we define the symmetric 
gradient by (B u)ij = diUj + djUi for i, j = 1,2,3; its deviatoric (trace-free) part is then 

= Du — - div ul, (1.3) 

O 

where I is the 3x3 identity matrix. The viscous stress tensor in each fluid is then given by 

§±(it±) := + p!± div u±I, (1.4) 

where p± is the shear viscosity and p!± is the bulk viscosity; we assume these satisfy the usual 
physical conditions 

p± > 0, p!± > 0. (1.5) 

The tensor P±(p±)I — §±(u±) is known as the stress tensor. The divergence of a symmetric 
tensor M is defined to be the vector with components (div M), = <9jM, r Note then that 

div (P±(p±)I — §±(u±)) = VP±(p±) - p±Au± - + // ± ) V div u±. ( 1 . 6 ) 

For each t > 0 we require that (p±,u±,r]±) satisfy the following equations: 


' d t p ± + div(p±u±) = 0 

p±(d t u± + u± ■ V«±) + VP±(/5±) - div§ ± (u ± ) = -gp±e 3 

dtjj± = u 3> ± - ui,±d yi r]± - u 2 ,±dy 2 r]± 

< (P + (p+)I - § + (ft+))n + = pat m n + - a + 'H + n + 

(P + (p + )I — S+(u+))n- = (P^(p-)I — §_(n_))n_ + 
u + = U- 
^u- = 0 


in £l±(t) 
in £l±(t) 
on r±(i) 
on T + (t) 
on T_(t) 
on r_(t) 
on 


(1.7) 


In the equations —gp±e 3 is the gravitational force with the constant g > 0 the acceleration of 
gravity and e 3 the vertical unit vector. The constant p a tm > 0 is the atmospheric pressure, 
and we take a± > 0 to be the constant coefficients of surface tension. In this paper, we let 
V* denote the horizontal gradient, div* denote the horizontal divergence and A* denote the 
horizontal Laplace operator. Then the upward-pointing unit normal of r±(i), n±, is given by 


n± 


(-V*r? ± , 1) 

Vr+Tv^p’ 


( 1 . 8 ) 


and TL±, twice the mean curvature of the surface r±(f), is given by the formula 


n± 


div* 


( v*v± \ 
\V 1 + \^ / *v±\ 2 ) ' 


(1.9) 


The third equation in (11.711 is called the kinematic boundary condition since it implies that the 
free surfaces are advected with the fluids. The boundary equations in (11,71) involving the stress 
tensor are called the dynamic boundary conditions. Notice that on T_(f), the continuity of 
velocity, u+ = ii—, means that it is the common value of u± that advects the interface. For a 
more physical description of the equations (11.71) and the boundary conditions in (11.71) . we refer 
to [20]. 

To complete the statement of the problem, we must specify the initial conditions. We assume 
that the initial surfaces r(0) are given by the graphs of the functions ? 7 ±( 0 ), which yield the 
open sets f2±(0) on which we specify the initial data for the density, /5±(0) : 12±(0) —>• M + , and 
the velocity, n±(0) : fl±(0) —> M 3 . We will assume that £ + 77 +(0) > r/-(0) > — b on T 2 , which 
means that at the initial time the boundaries do not intersect with each other. 
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1.2. Equilibria. We now seek a steady-state equilibrium solution to (11.71) with u± = 0, 77 - 1 - = 0, 
and the equilibrium domains given by the slabs 


= {y € T 2 x R | 0 < 2/3 < 1} and = {y € T 2 x R | — 6 < 2/3 < 0}. 


Then the system (11.71) reduces to the ODEs for the equilibrium densities p± = p±(y 3 ): 


d(P-fp-)) 

-- = ~9P~, 

a?/3 

-f+(P+(^)) = Patmi 

[p+(p+m = 


for 2/3 € (0,£), 
for 2/3 € (- 6 , 0 ), 


( 1 . 10 ) 


( 1 . 11 ) 


The system (11.111) admits a solution p± > 0 if and only if the equilibrium heights b,£ > 0, 
the pressure laws P±, and the atmospheric pressure p a tm fulfill a collection of admissibility 
conditions. These are enumerated in detail in our companion paper nu- For the sake of brevity 
we will not repeat them here, but we will assume that those admissibility conditions are satisfied 
so that an equilibrium exists. We remark that the equilibrium densities p are strictly positive 
and smooth when restricted to [— 6 , 0 ] and [ 0 ,K]. 

We denote the equilibrium density at the fluid interfaces by: 

Pi = P+(t), P + = P+( 0), P~ = P-(0). (1.12) 


Notice in particular that the equilibrium density can jump across the internal interface. The 
jump in the equilibrium density, which we denote by 

lp]'■= P+(0) ~ P-(0) = p + - p~, (1.13) 

is of fundamental importance in the analysis of solutions to (11.71) near the equilibrium. Since 
we are interested in the Rayleigh-Taylor instability, we assume [p] > 0, that is, the upper fluid 
is heavier than the lower fluid along the equilibrium interface. We refer to our companion paper 
nu for the analysis of the stable regime [[/?]] < 0. 

In studying perturbations of the equilibrium density it will be useful to employ the enthalpy 
functions. These are defined in terms of the pressure laws P± and the equilibrium density values 
via 

r z p' ( r ) rz p 1 ( r ) 

h + (z ) = / ——-dr and h-(z ) = / ——— -dr. (1.14) 

Jpi r J p- r 


1.3. Reformulation in flattened coordinates. The movement of the free surfaces T±(t) 
and the subsequent change of the domains £l±(t) create numerous mathematical difficulties. 
To circumvent these, we will switch to a coordinate system in which the boundaries and the 
domains stay fixed in time. In order to be consistent with our study of the nonlinear stability of 
the equilibrium state in HU , we will use the equilibrium domain as the fixed domain. We will not 
use a Lagrangian coordinate transformation, but rather utilize a special flattening coordinate 
transformation motivated by Beale [T] . 

To this end, we define the fixed domain 

fl = U with fl_|_ := {0 < ®3 < £} and := {—6 < X 3 < 0}, (1.15) 

for which we have written the coordinates as x € fh We shall write £_|_ := { 2:3 = £} for the 
upper boundary, £_ := { 2:3 = 0} for the internal interface and := { x 3 = —6} for the lower 
boundary. Throughout the paper we will write E = E_(_ U E_. We think of as a function on 
E-t according to rj + : (T 2 x {£}) x R + —»■ R and : (T 2 x {0}) x R + —>■ R, respectively. We 
will transform the free boundary problem in Q(t) to one in the fixed domain fl by using the 
unknown free surface functions rj±. For this we define 

f) + := V + rj + = Poisson extension of p + into T 2 x {^3 < £} (1.16) 


and 


:= = specialized Poisson extension of ?/_ into T 2 x R, 


(1.17) 
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where V± are defined in the appendix by (1A.4D and (1A.9I) . The Poisson extensions p± allow us 
to flatten the coordinate domains via the following special coordinate transformation: 

£l± 3 x ( xi,x 2 ,x 3 + b\fj + + b 2 rj_) := Q(x,t ) = (2/1,2/2,2/3) € S 2 ±(t), ( 1 . 18 ) 

where we have chosen b± = bi(x 3 ),b 2 = b 2 (x 3 ) to be two smooth functions in R that satisfy 

6i(0) = bi(-b) = 0, 61(f) = 1 and b 2 (l) = b 2 (-b) = 0 , b 2 (0) = 1 . ( 1 - 19 ) 

Note that 0(E + ,f) = r+(i), 0(E_,t) = T-(t) and 0(-,f) |e 6 = Id |s 6 - 

If 7) is sufficiently small (in an appropriate Sobolev space), then the mapping 0 is a diffeo- 
morphism. This allows us to transform the problem (11.71) to one in the fixed spatial domain 11 
for each t > 0. In order to write down the equations in the new coordinate system, we compute 


V 0 = 



and A := (V©” 1 ) 7 = 


1 

0 

—AK \ 


0 

1 

-BK 

(1.20) 

0 

0 

K ) 


: J' 

-1 


(1.21) 




(1.22) 


Here the components in the matrix are 

A = 0i0, B = d 2 0, J = 1 + d 3 9 , 1 

where we have written 

6 := bifj + + b 2 fj-. 

Notice that J = detV 0 is the Jacobian of the coordinate transformation. It is straightforward 
to check that, because of how we have defined fj- and 0 , the matrix A is regular across the 
interface 

We now define the density p± and the velocity u± on 11 ± by the compositions p±(x,t) = 
p±(Q±(x,t),t) an d u±(x,t) = u±(Q±(x,t),t)- Since the domains and the boundaries S-t 
are now fixed, we henceforth consolidate notation by writing / to refer to f± except when 
necessary to distinguish the two; when we write an equation for / we assume that the equation 
holds with the subscripts added on the domains £l± or To write the jump conditions on 
£_, for a quantity f = f±, we define the interfacial jump as 


If} f+\{x 3 = 0 } ~ I {3:3=0}' 

In the new coordinates, the PDEs f| 1 . 7 |> become the following system for ( p,u,r )): 


( 1 . 23 ) 


dtp — Kdt0d 3 p + div_ 4 (pu) = 0 in SI 

p(d t u - Kdfdd.iU + u ■ V^u) + SJ A P{p) - div^ § a(u) = -gpe 3 in SI 
dtp = u ■ M on E 

(• P(p)I - §a(u)W = PatmAl - (T+HM on E + (1-24) 

\P{p)I — §a(' u )1 N = v-HN on E_ 

[uj = 0 on E_ 

^U- = 0 on Efr- 

Here we have written the differential operators V_ 4 , div_ 4 , and A _4 with their actions given by 

(Va/)* := • Aijdjf , dlv A X := AijdjXt, and A A f := di v A X A f (1.25) 

for appropriate / and X. We have also written 

Jf ■■= (-dip, -d 2 p, 1) (1.26) 

for the non-unit normal vector to E(t), and we have written 


u)ij = AikdkUj + AjkdkUi , D^u = - - div^ ul, 

and S^4,±(rt) := p±H> A u + g!± div^4 ul. ( 1 . 27 ) 

Note that if we extend divyi to act on symmetric tensors in the natural way, then div_4 S A u = 
pA A u + (/r /3 + p')X A div A u. Recall that A is determined by 77 through (11.201) . This means 
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that all of the differential operators in (11.241) are connected to 77 , and hence to the geometry of 
the free surfaces. 


1.4. Perturbation equations. We will now rephrase the PDEs (11.241) in a perturbation for¬ 
mulation around the steady-state solution (p, 0,0). We define a special density perturbation 
by 

q = p-p-d 3 pd. (1.28) 

For the pressure term P(p) = P(p + q + d 3 p9), we expand it via the Taylor expansion: by (11.111) 
we have 


P(p + q + 8 3 p9) = P(p) + P'(p)(q + d 3 p0) +K = P(p) + P'{p)q - gpO + K, 


where the remainder term is given by 

CP+q+dspd 


n = 


I- 

J o 


(p + q + d 3 p0- z)P"{z) dz. 


(1.29) 


(1.30) 


Recalling (11.1311 . (11.191) . and (jl.221) . we find that 

- gp+0 = -P 19 V+ on s +> and l~gp9 1 = - {pj grj- on E_. (1.31) 

The equations (11.241) can be written as the following system when perturbed around the 
equilibrium (p, 0,0): 

dtq + divX(p + q + d 3 p0)u) — 8\pK9dt9 — Kdt6d 3 q = 0 in Ll 

(,P + q + d 3 pd)d t u + {p + q + d 3 p9)(—Kd t 0d 3 u + u ■ V^u) + pS7 a ( h'{p)q ) 

— div ^4 = —\7XP ~ g(g + d 3 p9)\7y\9 in Q 

dtp = u ■ AT on E 

(. P'(p)ql — §a( u ))N = PiggN — o + 'HN — 1Z + N on S + 

iP’{p)qI - S^(u)] M = [p] gpN + a-MM - [Tlj M on F_ 

[uj =0 on 

U- = 0 on Efe. 


(1.32) 


Remark 1.1. The special density perturbation q given by (|1.281) and the subsequent perturbation 
equations of the form (|1.32l) are crucial for our study in the stability regime in |llj . However, it 
is not essential for the instability regime in this paper. Indeed, we could consider p — p directly. 
We choose here to consider q in order to be consistent with the study in m- 


1.5. Main result. For a given jump value in the equilibrium density [p]| > 0, we define the 
critical surface tension value by 

°c ■= \P\g max{L^, L|}. (1.33) 

In our companion paper mi, we have proved the global existence of solutions decaying to the 
equilibrium state (0,0,0) in the problem (11.321) when cr_ > a c . The goal of this paper is to 
show that when cr_ < a c , the equilibrium state ( 0 , 0 , 0 ) is unstable in the compressible viscous 
surface-internal wave problem (11.321) . 

Our main result can be stated as follows: 


Theorem 1.2. Assume that [pj > 0 and a- < a c , where a c is defined by (11.331) . Let the triple 
norm |||-1| 00 be defined by (15.11) (with N > 3 an integer). There exist constants 9 q > 0 and C > 0 
such that for any sufficiently small 0 < t < 9 q there exist solutions {q b (t ), u L (t ), rf (t)) to (11.321) 
such that 

1 (^ (0)X(0),7/(0))1 00 < Cl, but ||XX)|| L2 > y. (1.34) 

Here the escape time T L > 0 is 

T L 1 log — 

A L 

where ^ < A < A with A the sharp linear growth rate defined by (12.341) . 


(1.35) 
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Remark 1.3. Theorem \1.2\ shows that the instability occurs in the L 2 norm of rThis high¬ 
lights the fact that the instability occurs at the internal interface. This also means that although 
our instability analysis works in a framework with some degree of regularity, the onset of insta¬ 
bility occurs at the lowest level of regularity. 

Remark 1.4. Our results can be readily applied to the compressible viscous internal wave prob¬ 
lem, i.e. the problem posed with a rigid top in place of the upper free surface. 

To our best knowledge, this work is the first rigorous result to address the nonlinear Rayleigh- 
Taylor instability for compressible viscous fluids with or without surface tension; for the com¬ 
pressible viscous internal wave problem, the linear instability was shown by Guo and Tice p 5 j. 
Theorem 11.21 together with our results in [TT| establish sharp nonlinear stability criteria for the 
equilibrium state in the compressible viscous surface-internal wave problem. We summarize 
these and the rates of decay to equilibrium in the following table. 



m <0 

[pi = o 

[p] > o 

a± = 0 

nonlinearly stable 
almost exponential decay 

locally well-posed 

nonlinearly unstable 

0 < a + 

0 < <7 < a c 

nonlinear ly stable 
exponential decay 

nonlinearly stable 
exponential decay 

nonlinearly unstable 

0 < <7+ 
a c = cm 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 

locally well-posed 

0 < <7+ 
a c < cm 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 


Note that our results do not cover the critical case: cr_ = a c . From [TO] we know that the 
problem is locally well-posed, but at the moment of writing it is not clear to us what the 
stability of the system should be. 

We mention some previous mathematical results concerning the Rayleigh-Taylor instability. 
For the inviscid Rayleigh-Taylor problem without surface tension, Ebin [3] proved the nonlinear 
ill-posedness of the problem for incompressible fluids, Guo and Tice [5] showed an analogous 
result for compressible fluids, and Hwang and Guo [3] proved the nonlinear instability of the 
incompressible problem with a continuous density distribution. For the viscous Rayleigh-Taylor 
problem, Priiss and Simonett m proved the nonlinear instability for incompressible fluids with 
surface tension in an LP setting by using Henry’s instability theorem, and Wang, Tice, and Kim 
mm established the sharp nonlinear instability criteria for the incompressible surface-internal 
wave problem with or without surface tension. 

Since linear instability can be established in the same way as that for the compressible 
viscous internal wave problem in [6] , the heart of the proof of Theorem 11.21 is the passage from 
linear instability to nonlinear instability. This is in general a delicate issue for PDEs since the 
spectrum of the linear part is fairly complicated and the unboundedness of the nonlinear part 
usually yields a loss in derivatives. Our proof is based on a variant of the bootstrap argument 
first developed by Guo and Strauss [3] . The main strategy of Guo-Strauss approach is to show 
that on the time scale of the instability, higher-regularity norms of the solution are actually 
bounded by the growth of low-regularity norms (in our case L 2 ). For our problem, the term 
[p] gr/_ along the interface is the cause of the instability; since it is of low order, we are led to 
use the Guo-Strauss bootstrap framework here. 

We encounter a number of mathematical difficulties in analyzing our complicated nonlinear 
problem, especially due to the fact that our problem is defined in a domain with a boundary. 
First, even to guarantee the existence of local-in-time solutions in our energy space, the initial 
data must satisfy certain nonlinear compatibility conditions that the growing modes to the 
linearized problem constructed in Section [2] would not satisfy. We employ an argument from 
Jang and Tice [9] that uses the linear growing mode to construct initial data for the nonlinear 
problem. 
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Second, because the spectrum of the solution operator for the linearized problem is compli¬ 
cated, we can only derive the largest growth rate for the linearized problem by using careful 
energy estimates as in [6]; this is in the context of strong solutions, which requires the initial 
data to satisfy the linear compatibility conditions. Such estimates would not be applicable to 
the nonlinear problem by directly employing Duhamel’s principle. To get around this issue, we 
provide the estimates for the growth in time of arbitrary solutions of the linear inhomogeneous 
equations in Section (3j clearly, the estimates can be applied directly to the nonlinear problem. 

The last difficulty is to derive the bootstrap energy estimates, a key step in showing the 
instability of r/_. We employ a variant of the energy method elaborated in our companion 
paper DU to derive these estimates in Section U There are new ingredients. First, we need 
to weaken the dissipation (see (14.211 for V % N ) due to the lack of the dissipation estimates of 
L 2 norm of q and r\. The missing terms are controlled by using the energy £% N instead. This 
is not possible in the global analysis of mi but is effective in the local-in-time framework of 
our instability analysis. Second, the term [pj gr /_ contributes negative energy in the unstable 
regime and thus the estimates of can not be obtained simultaneously with the other terms. 
We derive the estimates of r/_ by making use of the kinematic boundary condition, a transport 
equation for r/_. Our complete bootstrap estimate is recorded in Theorem 14.11 which shows 
that the stronger Sobolev norm II-Iqo of the solution is actually bounded by the lower-order 
norm ||r/_|| L2 . The bootstrap analysis allows us to finally prove Theorem 11.21 in Section [5j 

1.6. Definitions and terminology. We now mention some of the definitions, bits of notation, 
and conventions that we will use throughout the paper. 

Universal constants 

Throughout the paper we will refer to generic constants as “universal” if they depend on N, 
fi±, the various parameters of the problem (e.g. g , p±, p±, a±) and the functions p±. with the 
caveat that if the constant depends on a±, then it remains bounded above as either a± tend to 
0. For example this allows constants of the form gp + + 3 ct 2 _ + a + but forbids constants of the 
form 3 + 1/ a-. We make this choice in order to be able to handle together all the cases cr± > 0. 

We will employ the notation a < b to mean that a < Cb for a universal constant C > 0. 
Universal constants are allowed to change from line to line. When a constant depends on a 
quantity z we will write C = C(z) = C z to indicate this. To indicate some constants in some 
places so that they can be referred to later, we will denote them in particular by C\, C 2 , etc. 

Norms 

We write H k {Q±) with k > 0 and H S (Y>±) with s £ I for the usual Sobolev spaces. We will 
typically write H° = L 2 . If we write / € H k (U), the understanding is that / represents the 
pair f± defined on £l± respectively, and that f± € H k (Q±). We employ the same convention 
on E-i-. We will refer to the space 0 H l (Q) defined as follows: 

0 ={t)£ | [u] = 0 on E_ and V- = 0 on Ej,}. (1.36) 

To avoid notational clutter, we will avoid writing H k (£l) or H k (T,) in our norms and typically 
write only ||-|| fc , which we actually use to refer to sums 

\\f\\k = ll/+ll/7 fe (Q + ) + II/- lliT=(n_) or \\f\\k = II/+IIlU(E + ) + ll/-ll_f/ fe (E_) • (1-37) 

Since we will do this for functions defined on both and E, this presents some ambiguity. We 
avoid this by adopting two conventions. First, we assume that functions have natural spaces on 
which they “live.” For example, the functions u, p, q, and fj live on U, while g lives on E. As 
we proceed in our analysis, we will introduce various auxiliary functions; the spaces they live on 
will always be clear from the context. Second, whenever the norm of a function is computed on 
a space different from the one in which it lives, we will explicitly write the space. This typically 
arises when computing norms of traces onto E± of functions that live on U. 

Occasionally we will need to refer to the product of a norm of g and a constant that depends 
on ±. To denote this we will write 

7 INI l = 7+ b+llHUS+) + ^- Ik-ll■ ( L38 ) 

Derivatives 
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We write N = {0,1,2,...} for the collection of non-negative integers. When using space-time 
differential multi-indices, we will write N 1+m = {a = (ao,ai,. .. ,a m )} to emphasize that the 
0—index term is related to temporal derivatives. For just spatial derivatives we write N m . For 
a E N 1+m we write d a = d^d ^ 1 ■ ■ ■ . We define the parabolic counting of such multi-indices 

by writing |a| = 2ao + a\ + • • • + a m . We will write V*/ for the horizontal gradient of /, i.e. 
V*/ = d\fe\ + c? 2 /e 2 , while V/ will denote the usual full gradient. 

For a given norm ||-|| and an integer k > 0, we introduce the following notation for sums of 
spatial derivatives: 


v*7 


2 


E ii«“/ii 2 “ d 

ae N 2 

\a\<k 



E ii«"/ii 2 ■ 

aeN 3 

|a|<fc 


(1.39) 


The convention we adopt in this notation is that V* refers to only “horizontal” spatial deriva¬ 
tives, while V refers to full spatial derivatives. For space-time derivatives we add bars to our 
notation: 


v*/ 


2 


E ii«“/ii 2 and 

aGN 1+2 

\a\<k 


v fe / 


2 


E ii^iii 2 - 

aeN 1+3 

\a\<k 


(1.40) 


We allow for composition of derivatives in this counting scheme in a natural way; for example, 
we write 




£ ll<9“V*/|| 2 

aeN 2 

|a|<fc 


£ WfW 2 - 

aeN 2 

l<|a|<fc+l 


(1.41) 


2. Growing mode solution to the linearized equations 
In this section, we consider the linearization of (11,321) : 

' dtq + div(pu) = 0 
pdtu + pV (h'(p)q) — div§(ti) = 0 
dtJi = u 3 

(• P'ifi)ql - §M)e 3 = {P19V+ ~ o-+A*r/ + )e 3 
\P'{p)qI - S(u) 1 e 3 = ([p]] grj- + cr_A*p_)e 

H = 0 

^u- = 0 

We seek a growing mode solution to (12.11) of the following form: 

u(x,t) = w(x)e Xt , q(x,t ) = q(x)e Xt , rj(x',t) = ij(x , )e > 


in 

n 

in 

n 

on 

E 

on 

£+ 

on 

E_ 

on 

E_ 

on 


fj(x') 

)e xt 


( 2 . 1 ) 


( 2 . 2 ) 


for some A > 0 (the same in the upper and lower fluids), where x' = (xi,X 2 ). Substituting the 
ansatz (EH) into a, we find that 

q = — A -1 div(pu;) and fj = A _1 w; 3 |£. (2-3) 


By using (12.31) , we can eliminate q, fj from m and arrive at the following time-invariant system 
for w: 


X 2 pw — pV (h'(p) di v(pw)) — Adiv§(u;) = 0 in 12 

(-P'(p) di v(pw)I - A §(w))e 3 = (pigw 3 - a + Apw 3 )e 3 on E + 

< [-P'(p) di v(pw)I - A§(«;)] e 3 = ([p] gw 3 + cr„ A*w 3 )e 3 on E__ 

[u;]| = 0 on E_ 

W- = 0 on E;,. 


(2.4) 


Since the coefficients of the linear problem (12.41) only depend on the x 3 variable, we are free 
to make the further structural assumption that the x' dependence of w is given as a Fourier 
mode e lx for the spatial frequency £ = (£ 1 ,^ 2 ) E L~[ lr L x L^ 1 Z. Together with the growing 
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mode ansatz ( 12 . 21 ) . this constitutes a “normal mode” ansatz, which is standard in fluid stability 
analysis [2\. We define the new unknowns p,9,ip : {—b,£) —> R according to 

w\(x) = —ip{x 3 )e lx '6 W 2 (x) = —i 9 (x 3 )e lx '6 and w^(x) = ip{x 3 )e lx '6 (2.5) 


For each fixed 6 and for the new unknowns p{x 3 ), 0 {x 3 ),ip{x 3 ), and A, we obtain the following 
system of ODEs (here ' = d/dx 3 ): 


— (A pp')' + A 2 p + Ap |6 2 + (Ap' + A/r/3 + P'{p)p) J <p 
= ~6 [(Ap' + A/i/3) ip' + P'{p){pip)'} - 66 [A/i' + A/i/3 + P'{p)p\ 9 

— (A p9'y + A 2 p + A/i |6 + (Ap' + A/i/3 + P'{p)p) 9 

= -6 [(A/i' + A/i/3) ip' + P'{p){pipy] - 66 [A p! + A/i/3 + P'{p)p\ p 


- [( 4 A/i /3 + A/i') i/']' - p [/i'(p)(/n/0T + (^A 2 p + A/i |6 J 1 / 
= [(A/i' + A/i/3) (6<p + 6#)] + P [p'{p) ( 6 p + 60)] 
P+A(6^+ - p'+) = /i+A(6^+ - 0 +) = ° 


-(A/i' + A/i/3) (-0' + 6<P + 60) - -P'(p) ((pVO' + p(6<P + 6 9 )) 
-A/i (V>' - 6<p - 60) = (pip + C+ l^l 2 )^ 

M = M = M = [pA(6^ - p')l = [pA(6^ - 0')J = 0 
[(A/i' + A/i/3) (V’' + 6<p + 60)1 + [-P'(p) ((/#)' + p(6<p + 60))1 
+ [A/i (V'' - 6<p - 60)1 = -(IpI g-<r- 16 2 )V> 

^p = 9 = ip = o 


in {—b,£) 

in (—b,£) 

in {—b,£) 
at X 3 = £ 


( 2 . 6 ) 


at X 3 = l 
at X 3 = 0 

at X 3 = 0 
at X 3 = —b. 


We can reduce the complexity of the problem by removing the component 9. To this end, 
note that if p, 9, ip solve the equations (12.61) for £ € R 2 and A > 0, then for any rotation operator 
R G 50(2), (ip, 9) := R{p,9) solve the same equations for £ := R £ with ip, A unchanged. So, 
by choosing an appropriate rotation, we may assume without loss of generality that 6 = 0 and 
6 = 16 — 0- I n this setting 9 solves 


-(A pQ')' + (A 2 p + A/i |6 2 )0 = 0 

< 9{-b) = 9'{£) = 0 (2.7) 

M = IAp0'l = 0, 


which implies that 9 = 0 since we assume A > 0. Then the equations (12.61) are reduced to the 
equations for p, ip: 


— (A pp')' + A 2 p + Xp I6 2 + I6 2 (Ap' + Ap/3 + -P'(p)p) 

<P 


= - 16 [(Ap' + A/i/3 )iP' + P\p)(pipy] 


in (— b ,£) 

- [(4A/Z/3 + A/i') </>']' - P [h'{p){pipyy + (a 2 P + A/i I6 2 ) 

Ip 


= [(Ap' + Ap/3) 16 <p]' + P [P'{p) 16 <p] r 


in {—b,£) 

p+a(I6 - p'+) = 0 

-(Ap' + Ap/3 ){ip' + |6 p) ~ P'{p) {{pip)' + p 16 <P) 


at X3 = £ 

-Ap {ip 1 ~ 16 <P) = (PIP + ° + 


at X3 = £ 

M = M = [pA(l6 ip - <p')l = 0 

I(Ap' + Ap/3 ){ip' + 16 <p)l + lP'{p) {{pipy + P 16 <p)l 


at X3 = 0 

+ [Ap {ip' - 16 <p)l = -(IpI 9 - 0 - 16 2 )ip 


at X3 = 0 

p = iP = 0 


at X3 = —b. 


( 2 . 8 ) 


Solutions to (12.81) can be constructed in the same way as that for the compressible viscous 
internal wave problem in [ 6 ], so we will outline the procedure with minor modifications and 
refer some of the proofs to [ 6 |. 

It is not trivial at all to construct solutions by utilizing variational methods since A appears 
both linearly and quadratically. In order to circumvent this problem and restore the ability to 
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use variational methods, we artificially remove the linear dependence on A in (12.81) by introducing 
an arbitrary parameter s > 0. This results in a family (s > 0) of modified problems: 


- (W/ + A 2 p + sp |£| 2 + |£| 2 (sp 1 + sp/3 + P'(p)p) 

T 


= - l£l i(sp' + sp/ 3) if' + P'(p)(pif)'] 


in (—b,£) 

- [(4sp/3 + sp') y/f - p [h'(p)(p^yy + (a 2 p + sp lei 2 ) 



= [(sp' + sp/3)\Z\<p}' + plP'(p)\Z\v}' 


in (—b,£) 

A»+s(l£l^+ - T+) = 0 

-(sp 1 + sp/3)(if' + |£| ip) - P'(p) ((pif)' + yO ICI <p) 


at X3 = £ 

-sp (if - |£| v) = (pig + CT+ 


at X3 = £ 

M = M = IMIf 1 f - p')i = 0 

l(sp' + sp/3)(if + |£| 99)]] + iP'(p) ((pfi)' + p |£| 99)] 


at X3 = 0 

+ isp (if - |£| <p)J = -(IpI 9-0- \^\ 2 )f 


at X3 = 0 

(p = ip = 0 


at X3 = —b. 


A solution to the modified problem (12.91) with A = s corresponds to a solution to the original 
problem (12.81) . Note that for any fixed s > 0 and £, (12.91) is a standard eigenvalue problem 
for —A 2 , which has a natural variational structure that allows us to use variational methods 
to construct solutions. In order to understand A in a variational framework, we consider the 
energy functional 

E((p,i/j;s) = E 0 (ip,if) + sEi(<p,i/)) (2.10) 

with 

EoM) - — |{l ~y M V (o» 2 + "+ ^ |2 2 + PlS (m) 2 +1 fh'wmy+mv) 2 , e-n) 

= | J f* (w - l£l </0 2 + W - l£M 2 + !('/’' + Iflv) 2 ) + p'W'' + l£l ff- (2.12) 

and 

1 f* 

■%>V0 = 2 J pO 2 + V' 2 ), (2-13) 

which are both well-defined on the space ox $H l ((—b,£)) where 

0 H\(-b,£)) = {fi G H\(-b,£)) | <!>_ = 0 at x 3 = -b}. (2.14) 

Note that functions in this space automatically satisfy the condition [</>] = 0 at X 3 = 0. Consider 
the admissible set 

6 = € 0 H\(-b,£)) x 0 H\(-b,£)) | J(vA) = I}- (2.15) 

Notice that Eo((p,ip) is not positive definite for [p] > 0. 

The first proposition asserts that a minimizer of E in (12.101) over 6 exists and the minimizer 
solves (12.9p . 

Proposition 2.1. Let f and s > 0 be fixed. Then the following hold: 

(1) E achieves its infimum over &. 

(2) The minimizers are smooth when restricted to (—6,0) or (0,£) and solve the equations 
m with A 2 given by 

— A 2 = a(s) := inf E(<p,tp;s). (2.16) 

(¥M/>)ee 


Proof. A completion of the square and the fact that p solves (11.111) allow us to write 
h'{p) ((pip)' + p\f\p) 2 = P'(p)p(fi' + \£\<p) 2 -2gp\f,\if(p-2gptp'ip- gp'if 2 . 


(2.17) 
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We employ an integration by parts to see that 

/ -2 gpif'if - gp'if 2 = - [ g(pif 2 )'= -pig(if(£)) 2 + lp}g(if(0)) 2 . (2.18) 

J-b J-b 

We then obtain another expression for Eo(tp,if): 


eo(<p,i>) = + + 1 jy\p)p{^+ 1 ^^ ( 2 . 19 ) 

Notice that by further employing the identity —2 ab = (a — b) 2 — (a 2 + b 2 ) and the constraint on 
Jip), we see from (12.19|) that 

E(<p,if;s) > E 0 (p,if) > -2g\£\ [ pifp>-g |£| (2.20) 

J-b 

for any (ip,if) G 6. This shows that E is bounded below on 6. The results thus follow from 
standard compactness arguments, the variational principle for Euler-Langrange equations, and 
a bootstrap argument for smoothness. For more details, we refer to Propositions 3.1 and 3.2 of 
0 - □ 
In order to construct the growing mode solution to the original problem (12.11) we first need to 
ensure the negativity of the infimum (|2.16D . For er_ > a c , we always have that cr_ |£| — Jp] g > 0 
for any nonzero frequency £ G Lf lr L x Lf 1 Z, which implies E(p,if-,s) > 0 because of (12.111) . 
This then means that a(s) > 0, which suggests that no growing mode solution to (12.11) can be 
constructed when cr_ > a c , and in turn indicates that the system is linearly stable. In fact, 
in HU we have established the nonlinear stability of the compressible viscous surface-internal 
wave problem for the case <r_ > a c . However, when 0 < cr_ < cr c , for 0 < |£| < y/\p\g/aZ (it 
is interpreted that when cr_ = 0 this means 0 < |£| < oo), cr_ |£| 2 — [p] g < 0, and then it is 
possible for E(p,if; s ) to be negative. We denote this critical frequency by |£| c : 


leic : = 



( 2 . 21 ) 


Lemma 2.2. Let 0< |£| < |£| c . Then there exists sq > 0 depending on (j±,g, p,P,b,£, p,±, fj/±,\£\ 
such that for 0 < s < sq it holds that a(s) < 0. 


Proof. Since E and J have the same homogeneity, we may reduce to constructing any pair 
(ip, if) G o H l ((—b,£)) x o H l ((—b,£)) such that E(ip,if\s) < 0. We will take if with if(£) = 0 
and p = —if 1 / |£|, and we then further reduce to constructing any if G Hq((— 6, £)) such that 

E(if ; s) := E(-if'/ |£| ,if',s)= a -^ 2 ~ M 9 (if(Q )) 2 + \ f P'(PW + sE.i-if 1 / |£|, if) < 0. 

( 2 . 22 ) 

For a > 5 we define the test function tf a G -Hq((— 6, £)) according to 


lfa(x3) 



a/2 

a/2 


x 3 G [0, £) 
x 3 G (-6,0). 


(2.23) 


Simple calculations then show that 



\fT(b + £)T(a + 1) , . 

2T(a + 3/2) >: 


(2.24) 


where o a ( 1) is a quantity that vanishes as a —>• oo. We thus find that 

E(ifa, s ) = ^ |g| ~ M 9 + o a ( 1) + sC (2.25) 

for the constant C = E\(—if' a / |£| ,if a ), which depends on a, p,b,£, p±, p±, |£|. Since <7_ |£| 2 — 
Ip] g < 0, we may then fix a sufficiently large so that the first two terms sum to something 
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strictly negative. Then there exists so > 0 depending on cr±, g, p, P,b,£, n±, /j,±, |£| so that for 
s < so it holds that E(ip a -, s ) < 0. Thus a(s) < 0 for s < so- □ 

Remark 2.3. For a minimizer (p, ip) € & we have 

— ^ (^( 0)) 2 < a(s) < 0 , (2.26) 

which in particular requires that ip( 0) 0 and |£| 2 < g\p\/a-. 

As in Proposition 3.6 of [B], one can prove that a = a(s) is continuous and strictly increasing 
and that there exists s* € ( 0 , oo) so that 

5 = a _ 1 ((—oo, 0)) = (0, s*). (2.27) 

Arguing as in Theorem 3.8 of [ 6 ], we deduce the following. 

Lemma 2.4. For each fixed 0 < |£| < |£| c there exists a unique s € S so that A(|£|,s) = 
y/—a(s) > 0 and 

s = A(|£|,s). (2.28) 

Hence, we may now think of s = s(|£|) and we may also write A = A(|£|) from now on. In 
conclusion, we now have the existence of solutions to the system (12.61) . 

Proposition 2.5. For f € M 2 so that 0 < |£| < |£| c there exists a solution ip = tp(£,x 3 ), 
9 = 0(£,x 3 ), ip = ip(£,x 3 ), and A = A(|£|) > 0 to (12.61) so that ip(£, 0) 0. The solutions are 
smooth when restricted to (— 6 , 0 ) or (0,£), and they are equivariant in £ in the sense that if 
R G 50(2) is a rotation operator, then 

(p(R£,x 3 )\ (Rn Rvi 0\ (<p(£,x 3 )\ 

0(R£,x a) = U 21 R 22 0 9(£,x 3 ) . (2.29) 

\ip{R^,x 3 )J \ 0 0 1/ \ip(t;,x 3 )J 

Proof. We may hnd a rotation operator R S 50(2) so that Rf = (|£| ,0). For A = A(|£|) given 
in (12.28P, we define (<p(€,x 3 ),0(£, x 3 )) = i? _1 (^(|^| ,x 3 ),0) and ip(£;,x 3 ) = ip(\£\ ,x 3 ), where the 
functions </?(|£| , 2 : 3 ) and i/>(|£| , 2 : 3 ) are the minimizer of (12.161) . which solves the equations (12.91) . 
with s = A. This gives a solution to (12.61) . The equivariance in £ follows from the definition. □ 


To obtain a largest growth rate, we next show the boundedness of A(|£|). 
Proposition 2.6. For any 0 < |£| < |£| c , A(|£|) satisfies the bound 

bg [pj 


m\) < 


T- 


(2.30) 


Proof. For given |£| € (0,|£| c ), let (<p,ip) be the corresponding minimizer of E so that —A 2 = 
E(p, ip-, A). From (12.101) and (12.111) . we have E = Eq + A£i and Eq > — [pj g(ip( 0)) 2 /2. Hence, 


A E\ < — Eq < 


\9 


(mf 


On the other hand, since ip{—b) = 0, ip( 0) = J° b ip'dx 3 < Vb(J° b (ip') 2 dx 3 y/ 2 , and thus 

r0 ^ ^'-l^ , V>' + I ^ x2 

l-b 


(V>( 0)) 2 <b [ (ip') 2 dx 3 = b 
J-b 


dx 3 . 


2 2 

By further using the inequality: (A + B) 2 < 4(A 2 + B 2 / 3) for all A,BeR, we have 

(^(O )) 2 <b J ( ip ' - |£|v ?) 2 + ^(ip' + \f,\p) 2 dx 3 

= ~\ J b V- (itf ~ l^l ^) 2 + + |£|<// 2 ^) dx 3 < 


(2.31) 


(2.32) 


(2.33) 


Combining this with (12.311) . we deduce (12.301) . 


□ 
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Proposition 12.61 then allows us to define 

0 < A := sup A(|£|) < oo. (2.34) 

o<|£|<l£|c 

For cm > 0, only a finite number of spatial frequencies £ € (L^Z) x (Lf 1 Z) satisfy |£| < |£| c , 
so the the largest growth rate A must be achieved when 0 < cm < cr c . For cm = 0 it is not 
clear whether A is achieved. However, we can achieve a growth rate that is arbitrarily close to 
A, and so in particular A* is achieved, where 

0 < A/2 < A* < A. (2.35) 

We may now construct a growing mode solution to the linearized problem (12.11) . 

Theorem 2.7. Let A be defined by (12.341) and A* be defined by (12.351) . Then the following hold. 

(1) Let 0 < cm < a c . Then there is a growing mode solution to (12.11) so that 

IkWIU = e At \\ q m\k, |Kt)||* = e A *H0)|| fc , Ht)\\k = e At H0)\\k (2.36) 

for any k > 0. 

(2) Let cm — 0. Then there is a growing mode solution to (12.11) so that 

\m\\k = e A ik(0)|| fe , |Kt)||* = e A * t ||«(0)]| fc> Ht)\\k = (2-37) 

for any k > 0. 

Proof. Let |£| > 0 be so that A(|£|) = A for cm > 0 or A(|£|) = A* for cm = 0. Let (p = 3 ), 

6 = 0(£,x 3 ), if = 4>(f,x 3 ) be the solution to ()2.61) with A(|£|) as stated in Proposition 12.51 We 
then define q, u , and 77 according to (12.21) . (12.31) . and (12.51) . Then we have that q € H k (Q), 
u € 0 D and rj € H k (YL) for any k > 0 and ( q,u,rj ) solve the linearized problem 

m- Moreover, q,u,r] satisfy (|2.36l) or (12.371) . □ 


3. Growth of solutions to the linear inhomogeneous equations 

In this section, we will show that A defined by (12.341) is the sharp growth rate of arbitrary 
solutions to the linearized problem (12.11) . Since the spectrum of the linear operator is compli¬ 
cated, it is hard to obtain the largest growth rate of the solution operator in “L 2 —>• L 2 ” in the 
usual way. Instead, motivated by [6|, we can use careful energy estimates to show that e Ai is 
the sharp growth rate in a slightly weaker sense, say, for instance U H 2 —>• L 2 ”. However, this 
will be done for strong solutions to the problem, and it may be difficult to apply directly to 
the nonlinear problem due to the issue of compatibility conditions of the initial and boundary 
data since the problem is defined in a domain with boundary. We overcome this obstacle by 
proving the estimates for the growth in time of arbitrary solutions to the linear inhomogeneous 
equations: 


' dtq + divert) = G 1 

in Ll 


pdtu + pS7 (h'(p)q) — div§(ci) = G 2 

in Q 


dtp = u 3 + G 4 

on E 


( P\P)qI ~ §(u))e 3 = {pigp+ - cr + A* 77 + )e 3 + G\ 

on E + 

(3.1) 

\P\p)qI - S(u)j e 3 = ([p] gp- + cmA*? 7 _)e 3 - G 3 _ 

on E_ 


M = 0 

on E_ 


,u- = 0 

on E b , 



where G l, s are given functions. 

It will be convenient to work with a second-order formulation of the equations m- To arrive 
at this, we differentiate the second equation in time and eliminate the q and tj terms using the 
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other equations. The resulting equations for u read as 

pd tt u - pS7 {h'(p) div(ptt)) - di vS(d t u) = S 
(-P'(p) di v(pu)I - § (d t u))e 3 = {pigu 3 - CT + A*u 3 )e 3 + 0+ 
\-P'(p) di v(pu)I - S(d t -u)]] e 3 = ([pj gu 3 + cr_A*n 3 )e 3 - <3. 

[3t“l = 0 

dtu- = 0 


where 


and 


$ := -pV (h'ifiG 1 ) +8 t G\ 

©+ := -P' + {p + )G 1 e 3 - d t G\ + ( Pl gG% - u+A*G 4 )e 3 . 


in 

n 


on 

s+ 


on 

£_ 

(3.2) 

on 

£_ 


on 





(3.3) 



(3.4) 

)e 3 . 


(3.5) 


K(P+)G 1 ] e 3 - d t G 3 _ + 1 ~ ' 

Our first result gives an energy and its evolution equation for solutions to the second-order 
problem (13.21) . 

Lemma 3.1. Let u solve (13.21) . Then 
d / r - V'-' 


dt 


| \dtu\ 2 + |div(pu)| 2 + | |V*u 3 | 2 + / \u 3 \ 2 + 


/ an. 

4_ 2 


aM 


M' 


+ /£ 


Ddtu + Ddtu T — -(div dtu)I 


+ [ p'\divd t u\ 2 = I S • d t u - ( <5 ■ d t u. (3.6) 
J o J q, Je 


Proof. We multiply the first equation of (13.21) by dtu and then integrate by parts over fi. By 
using the boundary conditions in (13.21) . we obtain (13.61) . □ 

The variational characterization of A, which was given by (12.341) . gives rise to the next result. 

Lemma 3.2. Let u € o H l {Tl) n H 2 (Ll). Then we have the inequality 


I div (pu) | 2 + | |V*u 3 | 2 + yy |n 3 | 2 + 


■#M 2 


> a2 f -I ,2 A 

> —— / p \u\ 


P 


Du + Du t — - (div u)I 
o 


2 Jn 2 J n 2 

Proof. We take the horizontal Fourier transform to see that 


471-2 (X l div ^ pu ^ 2 + X f l v *' U3 l 2 + ^ 2 + 


+ p 1 |divn| 2 

5 M 


(3.7) 


W 


E 

^ei/^ZxLj 
+ y-l£| 2 + PiP 


ri h’ip) 2 

j ^ \ifipu i + + <9 3 (pu 3 )| 2 ete 3 


I^WI 2 + 


^-|e| 2 -[plp 


2 1 owl 2 
For f = 0, the term in the sum of (13.8|) is 


| 4 3 (<)| 2 - «M 


l%(0)| 2 

/l'(p) 


1^-3(0)| 2 + J ^~Y~ \d 3 {pv 3 )\ 2 dx 3 . 


i-b 2 

We expand the derivative term in the integral and integrate by parts to get 

r£ h’t p), a/= ., _ [* h'(p) 

i-b 2 

1 


(3.8) 


(3.9) 


J \d 3 (pu 3 )\ 2 dx 3 = j \d 3 pu 3 + pd 3 u 3 \ 2 dx 3 

1 


r - - - 

J 2 P '(p)p |d 3 ii 3 |^ - gpd 3 u 3 u 3 - -gd 3 p \u 3 \ z dx 3 

f b \p\p)p \d 3 h 3 \ 2 - P -f \um 2 + 1^3 (0)| 2 


(3-10) 
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This implies 

M| 4 3 (< )|2_9M| a 3(o)| 2 + J[ h ^\a 3{ pu 3) ^ X3 = \J[nP)mu 3 no. (3.H) 

Consider now the sum of (13.81) for fixed £ / 0, writing (p(x 3 ) = iui(f,,x 3 ), 9{x 3 ) = iu 2 {^,x 3 ), 
ip(x 3 ) = u 3 (t;,X 3 ). That is, define 

z(<p,9,ikZ) = a+ ^ 2 + Pl9 m)f + a ~^ 2 ~ Mg W0)) 2 

+ J + &p0 + {p^)'\ 2 dx 3 , (3.12) 

where ' = d 3 . The expression for Z is invariant under simultaneous rotations of £ and (ip, 9 ), so 
without loss of generality we may assume that £ = (|£| ,0) with |£| > 0 and 0 = 0. If <r_ > 0 
then we assume that |£| < |£| c as well. Then, using (12.101) with s = A(|£|), we may rewrite 

z (<P, 9 ,ip',€) =E((p,tp-, A(|£|)) - J // + |£| < p \ 2 

+ ( 3 ‘ 13 ) 

and hence 

z {< p , 6 ,^) >-y J p(M 2 + M 2 )- y/^'1'0'+ ICI^| 2 

/ / (k / -|^| 2 + ^-|^| 2 ^K + l^| 2 )- ( 3 - u ) 

Here in the inequality above we have used the following variational characterization for A, which 
follows directly from the definitions (12.16|) and (12.34p . 

EM; A(|£|)) > -A(|£|) 2 J(^,^) > -y J* p(\<p\ 2 + IVf)- (3-15) 


For |£| > £ c the expression for Z is non-negative, so the inequality (13.141) holds trivially, and so 
we deduce that it holds for all |£| > 0. 

Translating the inequality (13.141) back to the original notation for fixed £, we find 

+ °'- |€|2 T Wg |a s (0)| 2 + 


-TT- + i&pu 2 + (hlftU'i )] 2 dxs 


A 2 f e 


4 / /K.s. + ieA + W + f 


B 


(3.16) 


where 

B = Hu + Zlu T — -(div it)/. (3-17) 

O 

Taking sum of each side of this inequality over all 0 7 ^ £ g Lf 1 Z x Lf 1 ^, together with (13.lip , 
then proves the result. □ 


Now we can prove our main result of this section. We write 

Sg ■= IlG 1 !! 2 + \\d t G 2 \\l + ||5 t G 3 ||o + ||G 4 ||2 . (3.18) 

Theorem 3.3. Let ( q,u,r 7 ) solve (13.11) . Then we have the following estimates for t > 0: 

ll«(()li;<e 2A ‘(IMO)lll + l|9.«(0)llo)+ /Vf-'v^ll^Wlli*. (3.19) 

■Jo 

hMllo < V'(||«(0)|| 2 + ||9,«(0)|| 0 + h(0)|| 0 ) + [‘ ||G 4 ( S )|| 0 * 

Jo 


(3.20) 
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+ J J j e 2A ( s ■ '\V.V)| -2, „iri (/(/.,, 

llsWIlo £ e A '(ll“(0)ll 2 + l|S«»(0)ll 0 + ll«(0)||p) + f IlG'MH,* (3.21) 

Jo 

J J J e 2A ( s - r ) v 7 £g(t) ||<9 t u(T)H 1 drds , 


+ 


Proof. Integrating the result of Lemma 13.11 in time from 0 to t, and then applying Lemma I. '1.2 
we find that 


/>“' 2 +.U! 


Ddtu + Ddtu r — -(div dtu)I 
3 


Jo Jn 

< K 0 + f [ $-d t u- f [ <5-d t u 

Jo Jn Jo J e 

' h'(p) IJ; 


+ / /j,'\ div d t u\‘ 

Jn 


I div {pu)\ 2 J" |V*m 3 | 2 + ^y-|u 3 | 2 + 


In 2 

<K 0 +f [ $ ■ d t u - [ &-d t u 

Jo Jn 





0 J E 


A 2 


H—— / p\ u \ + 


A f P 

2 In 2 


Du + Z?rt T — -(div it)/ 


+ // |div u| 2 , 


(3.22) 


where 


R 0 = j f \dtu(0)\ 2 + I div (pu(0))| 2 


cr 


+ / — |V*u 3 (0)| 2 + / —|^3(0)1^ + 


g\pl 


l«s(o)|' 


(3.23) 


'E_ 


' E_ 


For notational simplicity we introduce the norms 


Ml l '■= [ p \ u \ 2 and IMI** '■= I ^ 

Jn Jn A 


Du + Du t — - (div u)I 


+// |divu| 2 (3.24) 


and the corresponding inner-products given by (■,■)* and (•,■)**> respectively. We may then 
compactly rewrite the previous inequality as 

\ WdMOWl + J o l|di u ( s )ll** ds <Ko + y 

where we have written 


< R o + — IKOII* + j IMOIL + fl(*) ( 3 - 25 ) 

t p pt 


fi(t) = $ -d t u- /©• d t u. 

Jo Jn Jo Jt, 

Integrating in time and using Cauchy’s inequality, we may bound 

rt 

2 


(3.26) 


A IM*)II** = A IK0)||** + A / 2 (u(s),d t u(s)) ick da 


io 


< 


A IM0)||** + [ ds + A 2 [ |Ks)|| 2 * ds 

Jo Jo 


On the other hand 


M ||«(t)||* = 2A(u(t),d t u{t)) ir < \\d t u{t)\\l + A 2 ||u(t)|| 2 . 

We may combine these two inequalities with (13.251) to derive the differential inequality 

dt IK*)II* + IMt)|| 2 * < K\ + 2A ^||«(t)||* + ||«(a)|lL ds ^j + jHO 


(3.27) 


(3.28) 


(3.29) 
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for K\ = 2Kq/A + 2 ||w(0)|| 2 . An application of Gronwall’s theorem then shows that 

ll«(i)lli + J‘ IkWIlL < e 2A ‘ ll«(0)lli + Ai(e 2A( - 1) + J f‘ (3.30) 

Now plugging (13.301) and (13.271) into (13.251) . we find that 

j \\dtu(t)\\l + ||ti(i)||^* < + A\\u{t)\\l + 2A ||u(a)||^ ds + jSj(t) 

<e 2M (2A\\u(0)\\l + K 1 ) + ^Sj(t) + / ^ f e 2A(t ~ s ^(s)ds. (3.31) 

A Jo 

Notice that 

A f\ 2A ^- s) Jj(s)ds = (e 2A(t " s) ) Sj{s)ds 

= --^-fj(t) + e 2A< i?(0) + J e 2At ' t ~ s ' ) d t S){s)ds (3.32) 


and 


£(°) = 0 and d t Sj = 3 


dtu — J 0 • d t u. 


(3.33) 


We then have, 

IK*)llL < e 2At (2A ||n(0)|| 2 + K\) + jJ^ e 2A(i ~ s) ■ d t u(s) - J ^ 0(a) • <9tu(s)^ ds. 

(3.34) 

By the trace theorem, 

Ki <lk(0)|| 2 + ||^n(0)|| 2 . (3.35) 


On the other hand, 


Iloilo + ll®llo ~ ll^li + ||5 t G 2 || 0 + ||i9tG 3 || 0 + ||G 4 || 2 — Eg- (3.36) 

So by Korn’s inequality (Proposition [AO]) and the trace theorem, ([3.341) implies (13.191) . 

Next we use the kinematic boundary condition dtT] = + G 4 and the trace theorem to 

estimate 

\\dtv(t)\\o ^ ll' u 3||//0( S ) + ||G 4 || 0 < 11 U 3 11 1 + ||G 4 || 0 . (3.37) 

This and (13.191) allow us to estimate 

\\v(t)\\ 0 < lb(0)|| 0 + [ \\dtv{s)\\ 0 ds < \\rj(0)\\ 0 + f \\G 4 (s)\\ ds + f ||« 3 (s)|| 1 da 
Jo Jo Jo 

<lli(0)ll„+ /'||G 1 ( S )|| 0 *+ /V(|K0)|| 2 + ||$«(o)|| 0 ) 

Jo Jo 

+ J J J e 2A ( s_r )-\/ Eg(t) ||5t'u(r)|| 1 drds, (3.38) 

which implies (13.201) . 

Similarly, we use the continuity equation dtq = — di v(pu) + G 1 to estimate 

II^Ho < ||div(ptt)|| 0 + ||G' 1 1| 0 < ||u||i + ||G' 1 || 0 • (3.39) 


We then deduce (|3.21|) as that for (13.201) . 


□ 
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4. Nonlinear energy estimates 


This section, the most technical part of the paper, is devoted to the nonlinear energy estimates 
for the system (11.32p . Our analysis here is similar to that found in our companion paper on the 
stable regime HU. The primary difference is that we will use slightly different versions of the 
energy and dissipation functionals in order to handle the fact that the internal interface makes 
a negative contribution to the energy and dissipation. 

For any integer N > 3, we define the energy as 


2N 

3=0 


4N-2j 


2N 

+ imiLv + Y ^ 1 

3 = 1 
2N 

+ I|V*t?|| In + \\v\\l N + Y |&tV 

3 = 1 


4AT-2j+l 


4JV-2j+3/2 


(4.1) 


and the “dissipation” as 


2N 


T> 2 N-=Y 


27V+1 


4JV-2/+1 


+ ll^gll^jv-i + Y, \ ^t q 


4N—2j+2 


3=0 j=2 

+ V 2 l|V*77|| 4Ar+1/ / 2 + V 2 11 ^7 11 47V -(-1/2 + H^ 7 ?ll4TV—1/2 + S IN 7 ? 


4iV+l 

E 

3 =2 


41V—2j+5/2 


(4.2) 


We also define 

3~2N '■= ll 7 ?ll41V+l/2 • (4-3) 

The surface tension coefficients a± are included in the definitions (14. ip and (14. 2p so that we will 
be able to treat the cases with and without surface tension together. It is noteworthy that the 
definition (14.21) of is different from the one introduced in {TT| for the nonlinear stability 
analysis: an< i ll r ?ll 4 iv-i /2 are n °t included in T>? iN here. This implies that to control 

or || 1 4 iv—l / 2 1 we have to use £ 2Ar ; mostly we will replace P 2JV in the estimates of some nonlinear 
terms derived in m by the sum + £® N . 

Our goal is to derive a priori estimates for solutions (q, u, rj) to (11.321) in our functional 
framework, i.e. for solutions satisfying £% N , Vf N , J+v < oo. Throughout the rest of this 
section we will assume that 

£%i v(?0 — — 1 (4-4) 


for some sufficiently small <5 > 0 and for all t € [0, T] where T > 0 is given. This assumption, in 
particular, will guarantee that the geometric terms introduced in Section 11.31 are well-behaved 
(see Lemma IA.3D . We will implicitly allow <5 to be made smaller in each result, but we will 
reiterate the smallness of 5 in our main result. Here is the main result of this section. 


Theorem 4.1. If sup 0<t<T £f N (t) < (5 2 for sufficiently small 5, then the following holds. For 
any e > 0, there exists C £ > 0 such that 


^+v(?) + J~2n(P) + [ F^Nds <C £ £f N (0) + J~2 n{ 0) + C £ f yj£f N ifDf N + £f N + J-2n) ds 

Jo Jo 


+ £ / (& 2 N + J~ 2 n) ds + C E / 11 ?7— 11 o ds 


(4.5) 


for all t € [0, T]. 


Theorem HU will be established by a series of energy estimates, elliptic estimates, and com¬ 
parison results and its final proof will be given at the end of this section. We start with the 
time differentiated version of problem (11.321) . 
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4.1. Energy evolution for temporal derivatives in geometric form. We will employ the 
form of the equations ([1.321) primarily for estimating the temporal derivatives of the solutions. 
Applying the temporal differential operator (/? for j = 0,..., 2N to ([1.321) . we find that 


dt{d{q) + di v A (pd{u) = F lj 

(p + q + dzp6)d t {diu) + pV A (ti(p)d{q) - di \ A E> A {d{u) = F 2 ’- 7 
dt (d1 p) = d{ u ■ Af + F 4, - ? 

( P'{p)d}ql ~ § A (d{u))Af = Pig&tV+N ~ °+^*{dtV+)N + F 
P'(p)cP t qI — E> A (d{u) Af = [[p| gd{p-Af + a-A*(d{p-)Af — F^ J on E_ 


& t u =0 
diu- = 0 


where 


F 1 ’ 3 = dj-F 1 — 22 CjdfAikdkipdl e u t ), 

o <i<j 


in 

n 

in 

n 

on 

E 

on 

S+ 

on 

E_ 

on 

E_ 

on 

S 6 > 


(4.6) 


(4.7) 


F?' J = &jF 2 + 22 C j { pAikd k ( df Ai rn dj e d m Ui) + ndfAikd} £ d k (Ai m d m Ui) 
o <l<j 

+ (p/3 + p)Aikd k (dtAimdi^ e d m ui) + (p/3 + p')dfA ik d J t ^ e d k (Aimd m ui) 
-pd$A ik dk{ti(p)di~ £ q) - df(q + d 3 p9)d t (&j.~ £ u) } , (4.8) 

Fi ’l = &{Fl + + 22 C j {/W<9f(A fiA ik )dl~ e d k ui + p+df(A fiAi k )d{~ £ d k Ui 

o<e<j ( 4 . 9 ) 

+(p'+ - 2p, + /3)d$(AfiAi k )di~ e d k ui + d £ t Afid{~ £ (pigp + - P\p)q - cr + A*p + )} , 

-Ft- = -dt F f- + 22 C J {d e t(MA ik )di~ e lpd k uij + d £ t (AfiAi k )d{~ £ \pd k Ui\ 
o <£<j 

+df(AfiAi k )di~ £ [(p' - 2p/3)d k uij + d £ t Afid{~\lp\ gp- - [P'(p)<?] + u_A*p_)} , 

(4.10) 


for * = 1 , 2 ,3, and 

F 4j = 22 CjdfAf ■ &j.~ £ u. (4.11) 

0 <£<j 

In the above, F 1 , F 2 and F 3 are defined by 

F 1 = d 3 pK9d t 9 + Kd t 9d 3 q - div^((g + d 3 p0)u), (4-12) 

F 2 = — (p + q + d 3 p6)(-Kd t ed 3 u + u ■ V A u) - V A U - g{q + d 3 p9)V A 9 , (4.13) 

F 3 = -77AC - <x+ div*(((l + |V*p + | 2 ) -1/2 - l)V*p+)Af, (4.14) 

and 

- F 3 = — pi} Af + a- div*(((l + |V*p _| 2 )- 1/2 - l)V*p_)AA (4.15) 


We present the estimates of these nonlinear terms F 1 ’- 7 , F 2j , F 3j and F 4 ’- 7 in the following 
lemma. 


Lemma 4.2. For each 0 < j < 2 N, we have 

I 2 . II T-,0 n | I 2 , II n 1 II 2 


^ lj '!lo + Iloilo + Iloilo + Iloilo £ £ 2N {V5n + S° 2N ) 


(4.16) 


Proof. The estimate is restated from Lemma 3.8 of m- Note that the appearance of (£ 2 n) 2 7n 
the estimate is due to the lack of llgll^ and IMI47V-1/2 7n the definition (14.2[) of P 2 n i we use 
£ 2 n to control them here. □ 

We now estimate the energy evolution of the pure temporal derivatives. 
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Proposition 4.3. It holds that 


/, 2 2 2 .2 \ F ~ J \ 

??+(*) n + cr ^*div(t) ) + Y1 

j =o V 1 o o o »/ io j=0 

< ^2Ar(0) + V£%N (V2N + ^2n) d- s + j 


\d 3 U 


ds 


7—No ds - 


(4.17) 


Proof. We take the dot product of the second equation of (14.61) with Jdju and integrate by 
parts over the domain fi; using the other conditions in (14.6p and some easy geometric identities 
involving J,A, and J\T, as in Proposition 3.1 of [IT], we obtain the following energy identity: 


1 d 


2 dt 

+ 


—+ ( / (. p+q + d 3 pO)J 
n 

2 


diu 


+ h'(p)J dlq 


+ / Pi9 
J Ex 


div-i 


+ Jp 


V*d 3 t r] 




J 


Jn 


Kdi 


u 


+ pJ di \^d{u 


= \ [ d t (J(p + q + d 3 pd)) dju +h\p)d t J d{q + [ J(h'(p)d 3 t qF 1 ’ j + dP t u ■ F 2 ’ j ) 

2 Jn Jn 

+ / -dt u ■ F 3 ' 3 + f p 1 gd 3 t p + F^ 3 - f aA^(dip)F 4,3 + f [p] (jdjrj_J\f ■ dju. (4.18) 
Jt, Jt, + ./e ./ e _ 

First, we may argue as in Proposition 4.3 of [7], utilizing Lemma lA.31 to estimate 

L‘ 


D \diu 


+ gj div_4 d^u 


~ Jn 2 


D °d 3 u 


div^u -Cy/E^V. 5 n . (4.19) 


We then estimate the right hand side of (14.181) for 0 < j < 2N. For the first two terms, we may 
bound as usual ||c^ J|| 

L°° < and \\d t (J(p + q + d 3 p9))\\ 

L°° < to have 

2 


[ d t (J(p 
Jn 


+ q + d 3 p0)) 


dfu 


< 

r^j 


F 


2 N 


dju 


+ 


d{q 


+ h'(p)dtJ 

~ y/e&eiN 


din 


(4.20) 


dju 


By Lemma 14.21 we may bound the and F 2, J terms as 

[ J{h , (p)d 3 qF 1 ’J + d(u ■ F 2 ’J) < dlq ||i ?1 ’ j |L + 

Jn 0 

~ \/^2N\J^2N (P 2 .N + ^2n)- 

For the F 3, J and F 4j terms, by Lemma 14.21 and trace theory, we have 




(4.21) 


L 


—d{ u-F 3, J+ f pigd^rj 

J E+ 


< 

rs_/ 


< 

rs_/ 


< 


di% 


diu 


H 0(E) 

+ 


S+ 

F 3 ’ j L+( di V+ 


J crA it (d 3 t p)F 4 ’ j 


+ (7 


A *<^77 




dtV+ Q + 0 A*d{ri ^ \Jz%n (P 2 N + ^2.n) 


y/Vht + ^2N\j^2N iP2N + ^2n) ' 


(4.22) 


For the last term, by the trace theorem and Cauchy’s inequality, we have 


/£_ 


Mgdiv-Af- diu < d 3 t p _ 


dtu 


H °(£_) 


< c £ 

r^j c 


dip- 


+ £ 


di 


u 


(4.23) 


for any e > 0 . 
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Consequently, employing Korn’s inequality from Proposition IA.4I in (14.191) together with the 
estimates (I4.20D (I4.23|) . taking e sufficiently small, and integrating (14.181) from 0 to t, we deduce 
that 


+ cr 


V* 


fyvit) + [ 

0 J o 


d^u 


~ ^2N\ 


ds. 


diq(t) 0 + 0 + 9tv+(t) 

r(0) + [ y/ & 2 N ^ 2 N + £%n) ds + f d{g 
Jo Jo 

Now taking j = 0 in (I4.24p . we have 

lk(i)llo + IK*)llo + ll 7 ?+( t )llo + cr ll V * r ?(Ollo+ [ \\ u \\l ds 

Jo 


ds 


(4.24) 


~ ^27v(0) + [ y /£%n (VZn + n) ds + f \\r]-\\l ds. 
Jo Jo 


For j = 1,... ,21V, the kinematic boundary condition, trace theory and the estimates 
show that 


dtV- 


< 


M 


H 0(E) 


+ IIf 4J_1 II 2 < 

~ II 110 ~ 


ar 1 


u 


+ N ('E > 2N + £-2n)- 


(4.25) 


1051 ) 


(4.26) 


Plugging (14.261) into (14.241) . by using — D we obtain 


d{q(t) + d 3 t u{t) + ^’ 77 +(t) 


+ u 


V, 


dtdit) + f 

0 Jo 


dfu 


ds 


<S. 


2 N 


(0) + f yj& 2 N (P 2 .N + ds + f 
Jo Jo 


lu 


(4.27) 


ds. 


Hence, by chaining together (|4.25l) and (I4.27|) . we get (|4. 171) . 


□ 


We remark that the energy identity in Proposition 3.1 of HD is slightly different from (I4.18|) . 
Unlike in Proposition 3.1 of HD , we do not employ the kinematic boundary condition in treating 
the last term in (14.181) because [p] > 0 ; if we did this, it would involve a negative term, 
— [/?]] g 11 77 — 11 0 , in the energy. As a result, for cm < cr C) the energy becomes non-positive definite, 
which is the cause of the instability. 


4.2. Energy evolution for horizontal space-time derivatives in linear form. We now 

estimate the energy evolution of the mixed horizontal space-time derivatives. It turns out to be 
convenient to rewrite the system (|1.32l) in a linear form such that the coefficients get fixed and 
that the elliptic regularity is readily adapted in later sections. The PDEs (11.321) can be also 
rewritten for (q, u, g) as 


dtq + div(pu) = G 1 

in 


pdtu + pS7 (h'(p)q) — div§(it) = G 2 

in S 7 


dtp = u 3 + G 4 

on £ 


(■ P'{p)ql ~ §(u))e 3 = (pigg + - a + A^g + )e 3 + G\ 

on £ + 

(4.28) 

\P'{p)qI - S(u) 1 e 3 = ([p] grj- + a-A*g_)e 3 - G 3 _ 

on £_ 


N = 0 

on £_ 


u- = 0 

on E 6 , 



where we have written the function G 1 = G 1,1 + G 1,2 for 

G 1 ’ 1 = Kd t 9d 3 q - uiAi k d k q , (4.29) 

G 1,2 = dlpKQdtO - qAi k d k ui - Ai k d k (d 3 pdui ) - (A ik - Si k )d k (pui), (4.30) 

the vector G 2 for 

G\ =-(q + d 3 pO)d t Ui + (p + q + d 3 p6)(Kd t 0d 3 Ui - uiAi k d k Ui) 

+ ^'AlkdkvA.lmdm'U'i “ 1 “ $lk$lm)dlzm^i 
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(4.31) 


T (/i/3 + fl )Ai k d k Al rn d rn Ui T (/i/3 + /i )(Ai k Al m &ik&lm)dkm'U , l 
- p(A a - S u )di(h'(p)q) - A u diK - g(q + d 3 pO)Audi6 
for i = 1,2, 3, the vector G+ = G ^ 1 + a+G 3 ^ 2 for 

G 3 ’l =n + (Audiu k + AkidiUi){Nk - 5k 3 ) + p+(Au - 5u)diu 3 + p(A 3 i - 5 3 i)diUi 
+ (p+ — 2fi + /3)AikdkUi(J\fi — 5i 3 ) + (p! + — 2fi + /3)(Aik — 5i k )d k ui5i 3 
+ pigr] + {Mi - 5 i3 ) - KNi + P'(p)q(5 i3 - Mf) (4.32) 

and 

Gj 2 = -A^ + (M - S i3 ) - div*(((l + (V^+l 2 )" 1 / 2 - l)V*r/+)A^ (4.33) 

for * = 1,2, 3, and the vector G 3 = G 3 / 1 + <r_G ^ 2 for 

-G 3 ’ 4 =(A« I/i^iifcl + Az [//9jUjl)(J\4 - <5 fc3 ) + (A - fe) [pd^I - (A - 5 3 i) {pdiui ] 

+ A [(// - 2g/3)d k uij (Mi - S i3 ) + (Aik ~ 5 tk ) [(//' - 2/t/3)5 fc uj] 5 i3 

+ Ip] gr)- (Mi -S i3 )- pZj Mi + lP'(p)ql (S i3 - Mi) (4.34) 

and 

G 3 ’ 2 = A* v _(Mi - 5a) +div*(((l + |V*r/_| 2 r 1/2 - l)V*r/_)^ (4.35) 

for * = 1,2, 3, and the function G 4 for 

G 4 = —uidirj — u 2 d 2 g. (4.36) 

We now present the estimates of these nonlinear terms G 1 , G 2 , G 3 and G 4 . Recall the 
notation V for space-time derivatives in (11.4(11) . 

Lemma 4.4. It holds that 

||v 4 Ar - 2 G 1 ||J +1|v 47 V_ 2 g 2 ||q + ||v 47 V ~ 2 g 3 || 2 /2 + ||v 4 Ar " 1 G' 4 ||5 /2 < e$ N (£% n + f 2N ) , (4.37) 


and 


V 4 


l G l 


lo 


+ V 4 


z d t G l ’ L \\l + ||v 47 V g 1,2 ||q + ||v 47 V - 1 g 2 ||q 


+ ii v 47 v_i g 3 ii ^ /2 +iivr- i G 4 nj /2 +nvr- 2 

< £ 2Ar (P 2 N + &2N + F 2 n) ■ 


^g 4 ||J /2 + cr 2 ||v 47 V g 4 || 2 /2 


(4.38) 


Proof. The estimates are restated from Lemma 3.3 of HP- The reason for the appearance of 
(“G27V ) 2 is the same as in Lemma 14.21 □ 

Next we present some variants of these estimates involving integrals of certain products. First 
we consider products with derivatives of G 4 . 


Lemma 4.5. Let a € N 2 such that lal = 4 N. Then 


and 


L 

L 


d a gd a G 4 


\J &2N {'E > 2N + £ 2 n) + 


aA*d'V"G 4 


^/ ^ 2 N^2N^2N 4” A /^9JV £') AT-F\ 


J 2N C '2N J ~ 2IV- 


Proof. The estimates are restated from Lemma 3.5 of El- 
Next we consider products with derivatives of G 1,1 . 

Lemma 4.6. Let a 6 N 3 such that |a| = 41V. Then 

J^h'(p)d a qd a G 1 ’ 1 < yJv° 2N + £$ N y/£$ N (V% + £» n + P 2N ). 

Proof. The estimate is restated from Lemma 3.6 of DU- 

We also consider a similar estimate involving weights and derivatives of G 1,1 . 


(4.39) 

(4.40) 

□ 


(4.41) 

□ 
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Lemma 4.7. Let a € N 1+3 with |ct| < 4iV and «o < 2 A J — 1. Then 

l (l + «“('>'(?) 9 )S° (fc'CpIG? 1 - 1 ) | < + 4 v +^)- 

(4.42) 

Proof. The estimate is restated from Lemma 3.7 of m □ 

We now estimate the energy evolution of the mixed horizontal space-time derivatives. 

Proposition 4.8. It holds that 


E (ll^Mllo + l|8““(*)llo + l|8V(*)llS + » l|8“l(*)ll?) + / E 

aSN 1+2 J ° a€N 1+2 

|q!|< 4A^ |a|<4A^ 

a 0 <2N-l a 0 <2N-l 


x u\\ 2 ds 


<£■ 


2 TV 


( 0 ) + f y/£f N (Dm + ££n + -^at) ds + / ||? 7 _ 
JO JO 


I4A-1/2 


ds. 


(4.43) 


Proof. Since the boundaries of are flat we are free to apply time derivatives and horizontal 
derivatives to the equations (|4.28j) . Let a € N 1+2 such that «o < 2iV — 1 and |a| < 4 Ah We 
apply d a to (14.2811 and argue as in Proposition 14.31 to obtain the following energy identity: 


~[ / h'{p)\d a q\ l + p\d a u\ z + I Pl g\d a g + \ 2 P I a\V*d a ri\‘ 


2 dt 


+ 


/E+ 
12 


^cxq3 


[ ^ |B°d a u| 2 + //|divd a u| 

Jn 2 

= [ h'{p)d a qd a G 1 + d a u • d a G 2 + [ -d a u ■ d° 

Jn Jt. 

+ f Pl gd a p + d a G\- [ aA*(d a r))d a G 4 + [ {pj gd a V .d a u 3 . (4.44) 

J S-i- J S J S- 

We hrst estimate the G 2 ,G 3 , G 4 terms in the right hand side of (14.441) . We assume initially 
that |a| < 4:N — 1. Then by the estimates (14.3811 of Lemma 14.41 we have 

fj a u-d a G 2 < ||a“u || 0 ||0“G 2 || O < y/V^yj£f N {V° N + £f N + P 2N ). (4.45) 

Similarly, the estimates (14.381) of Lemma 14.41 and trace theory show that 


[ d a u ■ d a G 3 

is 


r < ‘> I 


and 


L 


<ll^||^ (s) ||5“G 3 || 0 <||a“u|| 1 ||9“G 3 |l0 
~ \Z^2N\/ ^2N (PlN + N + -^2 n) 

(\\d a r,\\ 0 + a \\A*d a r,\\ 0 )\\d a G 4 \ 


(4.46) 


Pl gd a g + d a G\ - J aAJd a ri)d a G 4 


< 

f'XJ 


V^2N + &2N\J^2N 2N + &2N + J~2n) ■ (4.47) 


Now we assume that \a\ = 4 N. We first estimate the G 2 ,G 3 terms. Since «o < 2 N — 1, d a 
involves at least two spatial derivatives, and so we may write a = j3 + (a — f3) for some /3 € N 2 
with |/31 = 1. We then integrate by parts and use the estimates (|4.38l) of Lemma [4.41 to see that 


[ d a u-d a G 2 

— 

Jn 



d a ~ 0 G 2 

< 

d a+0 u 


d a ~ 0 G 2 




0 



< 


ll^tiHi ||V 


AN—lri2 


G 2 < 
llo ~ 


V^2N \J ^2N (P 2 N + ^2N + F2n) ■ 




(4.48) 
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Arguing similarly and using trace theory, we also find that 


L 


d* G 3 

= 

f d a+p u-d a ~ p G 3 

< 

d a+p u 


d a ~^G 3 



J E 



H- 1 / 2 (E) 



1/2 


<11^11, ||V 


47V-1.^3 


^ 3 |ll/2 ~ {^ 2 TV + ^ 2 .TV + ^ 2 n ) ■ 


For the G 4 term, we split into two cases: ao > 1 and ao = 0. If ao > 1, then d a involves at 
least one temporal derivative, so ||5 Q ry|| 3 y 2 < 11<9^°11 4 tv_ 2 cv 0 + 3/2 — ^ 2 N- This together with the 
estimates (I4.38() of Lemma 14.41 implies 


L 


Pl gd a r] + d a G 4 + - I aA*(d a rj)d a G 


/^4 


s+ 

< 


< II^L + ^I 


13/2 


I d a G 41 


11/2 


V^2N + ^2N \]£-21 V {^2N + &2N + -^Tv) . 


(4.49) 


If ao = 0, we must resort to the special estimates (I4.39|) (14.4011 of Lemma 1431 to bound, with a 
use of Cauchy’s inequality, 


[ Pl gd a r, + d a G\- [ aA^d a pd a G 4 
is± Jt, 


V£%n (P 2 N + £%n + ^ 21 v) • (4.50) 


We now turn back to estimate the G 1 term, and we recall that G 1 = G 1,1 + G 1,2 . For the 
G 1 ' 2 part, it follows directly from the estimates (14.381) of Lemma 14.41 that 


[ ti(p)d a qd a G 1 ’ 2 
Jn 


\\d a q\\o ||9°G 1 > 2 || 0 < v (P°N + £% n + T 21 v). (4.51) 


Now for the G 1,1 term we must split to two cases: ao > 1 and ao = 0. If ao > 1, then by the 
estimates (I4.38|) of Lemma 14.41 we have 


[ h'(p)d a qd a G 1 ’ 1 
Jn 


< ||a“g || 0 ||<9 Q G 1,1 || 0 < y/V^y/q N (P° 2N + £° n + T 2N ) . (4.52) 

If ao = 0, we must resort to the special estimates (14.411) of Lemma I4~6l to bound 

\/T , 27V + £%n ^ 2 n {^2N T £ 2 n + J~2i v) • (4.53) 


[ h'(p)d a qd a G 1 ’ 1 
Jn 


< 

r-^i 


Finally, for the last term in (14.441) . by the trace theorem and Cauchy’s inequality, since 
ao < 2 N — 1 and |a| < 4 N, we have 


j ^ M gd a r)-d a u 3 < \\d a g-\\_ l/2 \\d a u 3 \\ Hl/ 2 ^ ) 

<G £ ||^V||^_ 2 qo _i/ 2 + 6 ||^u||? 


(4.54) 


for any e > 0 . 

In light of (14.451) — (14.541) . we may now integrate (14.441) from 0 to t, apply Korn’s inequality, 
choose e sufficiently small, and use Cauchy’s inequality to find that 

\\d a q(t)\\l + \\d a u(t)\\l + \\d a V+ (t)\\ 2 0 + a ||V*<9^(i)llo + f \\d a n\\l ds 

Jo 

r(0) + y/£ 2 N (PZn + £ 2 N + J" 2 n) ds + 11 <9^° ? ?— 1147V— 2 Q !0 — 1/2 ds. (4.55) 


<£■ 


Now for ao = 0, summing (14.551) over such a gives 


E (ii^wiio + n^wiio + \\d a p+mi+* iiv^wng) + [ e 

N 
) 

^2Tv(0) + f \J£ 2N (JD 2 N + £Jn + J~2n) ds + [ I 
Jo Jo 


*w || 2 ds 


|a|<4TV 

a 0 =0 


|a|<47V 

Q 0 =0 


147V—1/2 


(4.56) 
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For cxq = j with 1 < j < 2N — 1, the kinematic boundary condition, trace theory, and the 
estimates (14.381) imply that 

|2 - - 2 • , . 2 


\K°V-\ 


< 

dj 1 u :i 

2 

+ 

di^G 4 



H 4N-2j-l/2(Y_) 



< 

< 




+ 


di^G 4 


AN—2j—1/2 


4 V—2(j —1)—5/2 


(4.57) 


+ £%N (J^2N + £-2N + J~2n) ■ 


y 4V-2(j 1)-3q}-1 u 

Plugging (14.571) into (14.551) . we obtain 

E (iia“9(0iiS+n«“<*(<)iiS + ii«“>i+(<)iiS+ffiiv.8“^(t)iiS) + [ E 


*u \\1 ds 


|oi| <4iV 
&o=j 

<£ 


2 N 


( 0 )+ 

J o 


2 N + £2N + J~2n) ds + 


\a\<4N 
&o=j 

'< v 4»- 20 -l)-3^-, u 2 


ds 


(4.58) 


since y/£% N < 1. Consequently, chaining (14.58j) and (I4.56P together leads us to (14.431) . □ 

4.3. Energy evolution for r] in transport equation. Note that the energy estimates of ?]_ 
are still missing in the energy evolutions presented in Sections 14.1114.21 We thus need to derive 
the estimates for r/_, and this can be done by revisiting the kinematic boundary condition, 
which is a transport equation for r]\ 

d t r] + u ■ V* 7 ? = U 3 in E, 

where u ■ V*r/ = u\d\r) + U 202 ??. 

Proposition 4.9. For any e > 0, there exists a constant C £ > 0 such that 
2 N „ rt 


(4.59) 


2 dtVtt) 

3 =0 


AN—2j 


<£ 


2 N\ 


and 


r( 0 ) + f \J£% N {F> 2 n + £ 2 N + ^ 2 n) ds 
Jo 

+ e [ £% N ds + C £ f Hv^ulli ds (4.60) 

Jo Jo 

F2n(J) < •7 7 2Ar(0) + C f \/ ^N-^iNds + e f F 2 N ds + C £ f ||V * JV «|| 1 ds. (4.61) 
Jo Jo Jo 

Proof. We first prove the estimates (14.601) . Recall that we have written dtrj = U 3 + G 4 . Applying 
d a for ol G N 1+2 with | a | < 41V to this and then taking the inner product with 9“ 77 , we obtain 

-jL\\d a v \\ 2 0 = [ d a r]d a u 3 + [ d a r]d a G 4 . (4.62) 

2 dt Jy. Jt. 

For the G 4 term, if ao > 1, then d a involves at least one temporal derivative, so the estimates 
(14.38 j) of Lemma 14.41 imply 


[ d a r]d a G 4 

is 


< Iloilo 


I^G 4 || 0 < i V 2 N + Zm + Kn)■ (4.63) 


If ao = 0, we use the special estimates (I4.39P of Lemma S3] to estimate 

[ d a r 1 d a G 4 < (Vf N + £f N + F 2N ) ■ 

Jt. 

On the other hand, trace theory and Cauchy’s inequality allow us to bound 

J d Q rjd a i 


(4.64) 


*7llo ll^ 3 || H o (s) < v^v ll^lli < eSi \n + \\Vi N u\\; (4.65) 


for any e > 0. Then the estimate (14.601) follows. 
































































26 


JUHI JANG, IAN TICE, AND YANJIN WANG 


To prove (|4.61|> . we define the operator J = \/l — A*. We apply J AN + 1 / 2 to (14.591) . multiply 
the resulting equation by J iN+l ^ 2 ri, and then integrate over £; using the standard commutator 
estimate, Sobolev embeddings on X, and trace theory, we find that 

1 d fiN = -^^-V*|J 4iV+1/2 r?| 2 + ^ (j m+1/2 u z - [j m+1 / 2 ,u] J 4N+1 / 2 V 


2 dt 


\ l^d lUl + d 2 n 2 )\J 4N+1 / 2 V \ 2 + f (j AN+l ' 2 u 3 - [j AN+1 ' 2 , u] • V*ry) 


~ II v * u IIl oo (e) 


+ 


j AN+1/2 


rj-AN + 1 / 2 r] 

V 


+ 


0 


0 


*'/ll L°°( E) 
AN 


jAN+1/ 2 u 3 
jAN+i/ 2 rj 


L2(S) 


IIV *u I 


L°°(S) 


jAN-l/2^^1 


0 


£ IMIffS(E) IMlLv+1/2 + II^Hli (! + Wvh) 11 ^ 114iV+l /2 
~ V^Zn^N + \Z?2N || V 4 


74N 4 2 i- 


Then the estimate (14.61(1 follows by using Cauchy’s inequality. 


(4.66) 

□ 


4.4. The evolution of energies controlling d 3 q. The energy evolutions presented in Sections 
I4.1H4.3I are not enough to get the full energy estimates by applying the Stokes regularity as in 
the incompressible case HSJI2] since we have not controlled divu. Motivated by Matsumura 
and Nishida [2], to control divu we introduce the material derivative of q in our coordinates: 


(4.68) 


Q ■= d t q - Kd t 9d 3 q + UjA jk d k q = d t q - G 1 ' 1 = - di v(pu) + G 1 ’ 2 . (4.67) 

We may then derive the following from (14.28|) : 

d 3 Q + p<9 3 (divu) = 0 3 G 1,2 - div(<9 3 pu) - d 3 pd 3 u 3 , 

pd t u 3 + pd 3 (h'(p)q ) - pAu 3 - (p/3 + p')d 3 (divu) = G\. 

By eliminating d 33 u 3 from the equations (|4.68l) . we obtain 

4p/3 + p\ , a \ 4 /V3 + /i' 0 ^ 1,2 , 1^,2 , 4 m/3 + /x' q 

h'(p)p 2 H + d ' S l ^ P ^ ~ - J 2 - d3 ° + ~p° 3 + h’(p)p 2 ddh 

- d t u 3 - 4 ^/ 3 + P (di v(d 3 pu) + d 3 pd 3 u 3 ) + ^(dnu 3 + d 22 u 3 - d 31 u 3 - d 32 u 2 ). (4.69) 

p- p 

In the light of (|4.67l) . we can view (14.69[) as the evolution equation for d 3 q. 

We now present the energy evolution of d 3 q. 

Proposition 4.10. For 0 < j < 2N — 1 and 0 < k < 4 N — 2 j — l, we have 

2 


k'<k 


o 


E pi N - 2j - k '- 1 d k 3 +1 ^h'(p)q) 

2j - k '- ld k'+idj (h'(p)q) * + E II \7i N - 2j - k '- l d% +1 diQ 


+ / X 

-2j-k'-l 

< £ftr(0) + f 

dl +l u 

Jo 


+ [ \! £ 2 n (PZn + £ 
Jo 


k'<k 


ds 


AN—2j— 


._ 1 + l|v^|li + E|| v ^- fe '^ 


k'<k 


k '+1 


ds 


\-J~2N)ds. ( 4 . 70 ) 

Proof. We first fix 0 < j < 2 N — 1 and then take 0 < k < dN — 2j — 1 and 0 < k' < k. Let 
a € N 2 so that |a| < 4 N — 2j — 1 — k! . Applying d a d\ dj to (14.691) . multiplying the resulting 
equation by d a d^ +l d{(h'(p)q) + d a d^ +1 d{(h'(p)Q), and then integrating over f2, we obtain 


I+ 11 + 111 = IV, 


(4.71) 
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where 


I = ladl'd! {h'(p)Q)^ d a dl' +1 di(h'(p)q), 

[ B a d k 'ai ( 4/i /' 3 + /i V h fd'fui 

J n 


II = I d a d k 3 &> y d 3 (h'(p)Q) ) 0“5 3 fc,+1 ^ (h'(p)Q) , 


/// = 


iV = 


[ d a dl' +1 di(h'(p)q)\ [ &*dg +1 dtih'(p)q)&*d* +1 di {h'(p)Q) , 

Jn Jn 

j + a°dl' +1 ei (h'(p)g )} 

x d°a$'di { + 4 g| + I^±fa 3 h'(m - d m 

{ P 2 P h'{p)p 2 

~ 4/i /3 + 11 (div(d 3 pu) + d 3 pd 3 u 3 ) + d n u 3 + d 2 2 U 3 - d 31 u i - d 32 u 2 ) 
P 2 P 


(4.72) 

(4.73) 

(4.74) 


(4.75) 


We will now estimate 1,11,111,1V. First, using the Cauchy-Schwarz inequality, we may 
easily estimate 


IV < 


d“d^ +1 dl(h'(p)p) 


o + £ 

k"<k' 


d a d k " +1 diQ 


d^G 1,2 


4N—2j 


+ 


%G 2 


4N—2j—l 


+ 1 


k"<k' 


+ 


&i +l i 


d a dfd{u 


+ 


V*V<^u 


(4.76) 


+ V 

4A-21-1 

J k"<k' 

For the last term in III we recall the definition of Q from (14.671) in order to rewrite 

[ d a d k 3 ' +1 di(h r (p)q)d a dl' +1 di ( h\p)Q) 

Jn 

= [ d a d k 3 +1 &[(h\p)q)d a d k 3 +1 di (h'{p){d t q- G 1 ’ 1 )) 

Jn 

= \if |a“9i' +1 ^(ft'(p)9)( - / 9”8 3 ‘ ,+1 9|{/.'(p)«)9“9 3 t '+ 1 fl 2 ' (A'(p)G 1 - 1 ). (4.77) 

For II we estimate by expanding with the Leibniz rule: 

7 4/i/3 + p,' • 2 


77 > 7c ^'(/>)p 2 

For /, we have 

I > 


d a d 3 +1 d{ (h'(p)Q) -C d a d k ' +1 dt (h'(p)Q) I d°d 3 ' cf t Q 


k"<k' 


o 

(4.78) 


jf (h'(p)Q) (Tdl'+ l Si(h'(p) q ) 


-c 


d a d k 3 +1 &t(h\p)q) ]T I d a d k "d}Q 


k"<k' 


Id f dp/2, + p' 


d a d k ' +1 di(h\p)q ) 


2 df d'(p)p 2 
- f a A ^ p 2 d a d k 3 +1 di(h\p)q)d a d k 3 + 1 di {h'{p)G^) 


-c 


d a d k 3 +1 &>{h\p)q ) £ d a dfd J t Q 


k"<k' 


(4.79) 
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Combining the estimates (14.761) (14.791) with (14.711) . applying Cauchy’s inequality in order to 


absorb the term 


d a d k 3 +1 dl ( h!(p)q) 


we arrive at the inequality 


onto the left, and then integrating in time from 0 to t, 


d a dl' +1 di(h'(p) q )\ 2 + f d a d k 3 +1 dj{h\p)q) V d a d k 3 +l &> (h'(p)Q) 
1 0 Jo 0 

<£m(0)+ f E | d a ^'&lQ 

J 0 /„// ^ uf 

+ L I L 
l 


ds 


(4.80) 


k"<k' 

4/r/3 + p'\ a k t 


h'{p)p 


d a d k ' +1 di(ti(p)q)d a d k ' +1 di (h'(p)G 1 - 1 ) 


d(G 1,2 

2 

+ 

&{G 2 

2 

+ 

d} +1 u 

2 

+ 

v 47V-2 j-k'Qj u 


47V—2 j 


47V—2j —1 


1 

to 

^o. 

1 



+ 


Owing to the Leibniz rule and the properties of p, we may estimate 


fc'+i 


ds. 


d k ' +1 d a djQ 


< 


ti(p)d 3 +1 d a d{Q 


c^' +1 a“c^(ft'(p)e)|| + XI ||#>ate|L- ( 4 -»i) 

k"<k’ 

Combining this with (14.80j) and summing over all a. with |or| < 4 N — 2j — 1 — k', we deduce that 

+ f ' +1 di (h'(p)q) 2 + Vi N - 2j - k '- l d k ' +l d]Q 2 ds 

Jo 0 0 

<£m(0)+ f E I Vi N -^- k '- 1 d k "d{Q 2 ds 
Jo k"<k ’ 1 0 

+ f E [ ( 1 + \u-t 2 ' )d a d k 3 +i di(h'(-p)q)d a d k 3 +i di{h\p)G^) 

J ° |a|<47V—27 —1—fc' Jn V h {P > P J 


ds 


+ 


f 

Jo 


diew 


+ 


d}G 2 


4N—2j—l 


+ 


di +1 u 


47V—2j —1 


+ 


Vi N - 2j ~ k 'd}u 


k '+1 


ds (4.82) 


4N—2j 
for each 0 < k' < k. 

Finally, we will estimate the nonlinear terms in the right hand side of (14.821) . We use the 
estimates (14.38 j) of Lemma 14.41 to estimate, for 0 < j < 2N — 1, 


3'G 1 ’ 2 


+ 


4N—2j 

Then we use Lemma 14.71 to bound 
4/r/3 + p! 


diG 2 


47V—2j —1 


< £% n {T>2N + £2 TV + ^2n) ■ 


(4.83) 


/ 

Jn 


1 + 


h'(p)p 2 


d a d k +1 di(h'(p)q)d a dI +1 d{ {h'(p)G^) 

~ \/'£ > 2N + £"2N ^2N {P 2 N "t" ^2N J~2n') ■ (4.84) 


Plugging the nonlinear estimates (I4.83|) and (|4.84|) into (I4.82|) then yields that every 0 < k 1 < k, 

2 

o 


v 4IV-2 j - fe '- 1 ^' +1 ^ ( ^ ( - )g) 


+ 


v 4N-2j-k'-l d k'+lQj (^(^g)! + ||v^-2j-fc'- 1 a fc'+ 1 ^g 


ds 
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AN-2j-l 


+ 


v m-2j-k'oiu 


k '+1 


+ V £%n (^2 n + £%i v + J~ 2 n) ds (4.85) 


since ^/£ 2N — 1- We recall the notation Q in (14.671) . We may use the estimates (14.381) of Lemma 
14.41 to obtain the bound 

||' Vi N Q\\l < II vr div(pu)||o + || V^G 1 - 2 ^ < || V^v-Wl + £% N (V% n + £% n + F 2N ) . (4.86) 

Then a standard induction argument on (j4.85[) . together with (I4.86p . yield (14.701) . □ 

Note that the novelty of Proposition 14.101 is twofold. First, it presents the energy estimates 
of dsq. Second, it provides the dissipation estimates of dsQ and hence dsdiv(pu) by (I4.67p . 
These are crucial for improving the horizontal energy and dissipation estimates derived in the 
previous section into the full ones in later sections, respectively. 


4.5. Combined energy evolution estimates. Now we chain the results in Sections 14.1114.21 
and 14.41 with the elliptic regularity theory of a certain Stokes problem into an intermediate 
energy-dissipation estimate. 

We first derive the elliptic estimates. We deduce from (14.281) that 
div(pu) = G 1 ’ 2 - Q, 

- ^ Au - + M V div u + V (ti{p)q) = -G 2 - d t u. 

p p p 

Direct calculations give the form of the Stokes problem we shall use: 


-pA + V (. h'(p)q ) = ^G 2 - d t u- p 

| P/ 3 + p' 

P 

div (G 1,2 - Q - 2dspus) 

u = u 

We now prove the Stokes estimates. 


- P ( 2<9,3 Q J dsu + dss 
V (J (G 1 ’ 2 - Q-dspus) 


u 


in £l± 

in D-i- 
on <90-i-. 


(4.87) 


(4.88) 


Lemma 4.11. Fix 0 < j < 2 N — 1. Then for any 1 < k < 4 N — 2 j, 


v 47V-2 j- k d j u 


< 


fc+1 


+ 


Wt N -2 j - k di (h'(p)q) 


d{ +1 u 

2 

+ 


4IV—2j—1 


Vi N ~ 2 j ~ k diQ 


2 

fe-1 
V 74 Af „,|| 2 


+ || V* u\\ 1 + £f N ^ 2 N ^n) ■ (4.89) 


Proof. We first fix 0 < j < 2N — 1 and then take 1 < k < 4 N — 2 j. Let a € N 2 such that 
|ck| < 4 N — 2 j — k. We apply d a df[ to the equations (14.881) in D-t respectively; then the elliptic 
estimates of Lemma IA.6I with r = k' + 1 > 2 for any 1 < k! < k and trace theory allow us to 
obtain the bounds 


d a d J t u 


k '+1 


+ 


Vd a d} (h'(p)q) 


2 

< 

7 . 

d a di(-) 

2 

+ 

k'—l 

\pj 

k '+1 


Vd a di (h'(p)q) 


k'-l 


< 

r^j 


d a djG 2 


k'-l 


+ 


d a d{ +1 u 


+ 


d a d{ G 1,2 + d a diQ 

k' 


k'-l 


+ 


+ 


d a cf t u 


k> 


d a d{u 


< 


d a d J t G 1,2 


+ 


d a d{G 2 


k-1 


+ 


d Q di 


+i 


H k '+ 1/2(S) 
2 

k -1 


d Q d{ Q 


+ 


d a cf t u 


k' 


+ 


V k 'd a d{u 


H 1 / 2 (E) 


< 

t~KJ 


die 1 * 2 

2 

+ 

die 2 

2 

+ 

d{ +1 u 

2 

+ 


AN-2j 


AN—2j—1 


4AT-2J-1 


Vi N - 2j - k djQ 
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+ ll^* 4 NI? + 


d a d{u 


k> 


(4.90) 


A simple induction based on the above yields that 


d a d}u 


fe+i 


+ 


X7d a dlq 


k -1 


< 


djG 1 > 2 


4N-2j 


+ 


dlG 2 


+ 

di +1 u 

2 

+ 



4AT-2j-l 


v 4N-2 j-kQj Q 


4N—2j—l 

+ ||v^||' 


Finally, we use the estimates (14.381) of Lemma 14.41 to have 


diG 1 ' 2 


m-2j 


+ 


did 2 


< 


4JV-2j-l 


£%N {'D‘2 N + £%.N + ^2N) 


We then sum (14.911) over such |a| < 4IV — 2j — k to conclude (14.891) . 


(4.91) 

(4.92) 

□ 


We will now combine the energy evolution estimates of Sections l4.1H4.2l with the d^q estimate 
of Section 14.41 and the estimates of Lemma 14.111 The full dissipation estimates of u will be 
obtained, and also some estimates of q will be improved along the way. To do so, we first 
introduce some notation. We write 


and 


2N 

£■2N = llvf^llo + Iloilo + E |^+ 

1=0 

®2N = IIV^IZ 


2 

4N-2j 


2N 

E 

1=0 

2 

4N—2j 

(4.93) 



(4.94) 


for the various terms appearing in Propositions l4.3l and l4.81 Similarly, for integers 0 < j < 2N—1 
and 0 < k < 41V — 2j — 1 we write 


■= E I vi N -v- k '-'el'+'di(h'(p)g) 
k '=0 


(4.95) 


and 


<B& := E I Vi n -^- t '- I d^ I 9; (h'(p)q) o + E ||wi' v - 2 >-'‘ , 8|{fc'(p) 9 ) 

k '=0 k '=1 

+ ||vrQ||o+ E I 


2 

fc'-l 


fc'=0 


(4.96) 


where Q is defined in (14.671) . In addition, we introduce the following intermediate energies: 

2N 0 2N n 2N 

V74AL.II 2 


and 


£%n := || v ^|| 0 + E H 9 

l=o 

2 N 

V2N := E 1^ 
1=0 


4N-2j 


1=0 


+EII^+ 4Ar _ 2j +c7 Ell v ^ r ? 


1=0 


2 

4N—2j 


U 


71 — 1 


4N-2j+l 


+E 


1=0 


4N-2j-l 


(4.97) 


(4.98) 


The rest of the section is devoted to the derivation of the energy bounds for £% N and T> 2 n based 
on the evolution equations for ® 2 AL 2l 2W -, ® 2 )v 

Proposition 4.12. Let £% N and T> 2 N be as defined by (14.971) and (|4.98l) . Then we have 


p 

Cs‘ 


2 N 


(t) + f T> 2 Nds < f 2 7v(0) + f \J& 2 N (^ 2 A r + £- 2 N + ^ 2 N ) ds + f \\t]- / 2 ds. (4.99) 

J 0 Jo Jo 
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Proof. First, we sum the result of Proposition 14.31 with the result of Proposition 14.81 this yields 
the estimate 


e 


2 N 


(t) + f '£>2Nds < £f N ( 0 ) + f \J£f N (T>2 n + £%n + -Fi n) ds + f ||? 7 — llliv—1/2 ds, ( 4 . 100 ) 
Jo Jo Jo 


where <£? 1N and £>2 jv are as defined in (|4.93[) and (14.941) . 

Next, for 0 < j < 2N — 1 and 0 < k < 4N — 2j — 1, we may combine the results of Proposition 
14.101 and Lemma 14.111 (summed over 1 <k’<k) with (|4.86D to see that 


%2N + [ ~ ^2Af(0) + 

Jo 


dj +1 u 


2 

AN-2j- 




ds 


k '=1 


+ f ®2 N + \/£%N i^2N ^2N + J~2n) ds, 
Jo 


(4.101) 


where 21 J 2 ’^ and *8^ are given in (|4.95[) and (|4.96l) . Note that Lemma [4.111 used here is to 


control the term 


then we have 


In turn, we have 


v 4IV-2 j-k'gj 


u 


2 

k'+l 


in the right hand side of (14.701) . If we write 


■= ||V^Q||o+ Vi N ~ 2j - k - l cP t Q 


2 

fc+i 


&& < l|v^Q||o + E I 

k '=0 


— ■°2N- 


(4.102) 


(4.103) 


%2N + [ &2N^ S ~ v(0) + [ &t +1,1 

J 0 Jo 


4N—2j— 1 


E hi N - 2j ~ k 'diQ 


k'= 1 


ds 


+ / ®2 N + \/£%N 0^2N + £%N + J~2N ) ds. 
Jo 


(4.104) 


A standard induction argument on the above yields 
4A-2J — 1 1 4N—2j—l 

'0 


r 

d{ +1 u 

2 

+ 

Vi N - 2j djQ 

Jo 


4iV-2j-l 



E ^2A + [ E ^2N^ S ~ ^27v(°) + [ 

■ " k=0 J ° 

+ [ ® 2 iV + \/£%N i^ 2 N "P £fw + J~2n) ds 
Jo 

( 0 )+ / 

Jo 


ds 


dl +l u 


AN—2j—1 


+ + \/£%n 0^2N “P £&v + J~ 2 n) ds, (4.105) 


where we have used (14.861) to derive the second inequality. 


Note now, using the definition of £j^, that 


VrS||o+ %Q 


AN—2j—l 

< cJ’ k 


4N—2j 


E ^2 


/c=0 


(4.106) 


Using Lemma 14.111 with k = 41V — 2j, we then have that 


dfu 


+ 


4N-2j+l 
AN—2j—l 

~ E ^2A + 

k=0 


Vdi(ti(p)q) 


dj +1 


u 


4N—2j—\ 

+ £) 2 A + & 2 N i^ 2 N + £ 2 N + J~ 2 n) ■ 


AN—2j—l 


(4.107) 
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Hence (14.1051) implies 

41V-2J-1 f 

E a X + / 

r\ JO 


k=0 

<£■ 


dji 


2N 


(o) + r 

Jo 


di +1 u 


AN—2j+l 
2 


+ 


Vd{(h\p)q) 


AN—2j—l 


ds 


+ + \J£%n 0^2n d~ £%n + J~2n) ds , (4.108) 


4AT—2j —1 

for all 0 < j < 21V — 1. A standard induction argument on the above yields 
2AT—1 4AT— 2j-1 t 2N-1 


Y X] + j Y 

JO j=0 


j =0 fc=0 

ib C 2AT 


dju 


AN-2j+l 


+ 


VdJ t (h\p)q) 


AN—2j—A 


ds 


( 0 ) + f '&2N + \J£%N 0^2N + <?2A T + J~2N) ds. 
Jo 


(4.109) 


Consequently, a suitable linear combination of (14.1001) and (14.1091) gives 
2JV—1 AN-2j-l t 


^2n(0 + Y/ X/ ^2A + f £ > 2Nds 

j =0 fc=0 

<^ 2 Ar( 0 )+ f y 1 £% n {T> 2 N + £%N + ^ 2 n) ds + f ||r/_| 
Jo Jo 


AN-1/2 


ds. 


(4.110) 


Note that 


2N—1 AN—2j—l 2N-1 AN-2j-l k 

E E a l£ - E E E 

j= 0 /c=0 j=0 /c=0 /c'=0 

2AT—1 4N—2j — l 27V—1 

= E E |vi A '- 2 '-*- 1 %8?A'(p)9)| t - J] ||M(k'(?)g) 

j=0 fc=0 j=0 

Since 

d 3 dlq = -J— d 3 dl(h'(p)q) - d 3 {h\p))dlq 
and h'{p) is smooth on [—6,0] and [0,-f] and bounded below from zero, we may estimate 


2 

o 

2 

4JV—2j —1 


:=Z. (4.111) 


(4.112) 


d 3 d(q 

2 

< 

0 ~ 

d 3 d{ (h'(p)q) 

2 

+ 

0 

d{q 

2 

< z + 

0 ~ 

dlq 


o 


zz+m^qi 


0 


(4.113) 


and similarly 


d 3 d}q 

2 

< 

rsj 

l 

d 3 df(ti(p)q) 

2 

+ 

i 

d{q 

2 

< z + 

l 

VI dlq 

2 

+ 

0 

d 3 diq 


2 

i —1 


<Z+\\wi N - 1 q\\l+ d ^tQ ,_ x (4-114) 


for 1 = 1,..., 41V — 2j — 1. A standard induction argument then yields 

4JV—2j — 1 


d 3 d J t q 


AN—2j— 


l < Y dsd t q 


i =0 


<z + llvf^d 


On the other hand, 


dtq 


AN—2j 

so summing (14.1151) and (14.1161) yields 

2N 

Y 

3=0 


<||V^||o + 


d 3 cP t q 


AN—2j — l 


AN—2j 


<Z+\\Vj N q 


* ^ 110 


We then deduce (I4.99P from (14.1101) and (14.1171) . 


(4.115) 

(4.116) 

(4.117) 

□ 
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4.6. Full energy estimates. In this section, we will derive our ultimate energy estimates. 
First, we combine the results in Sections 14.31 and m We define 

2 N 


2 

4.N-2 j 


(4.118) 


£■2 Tv(^) + J~2N 


(t) + f 

JO 


&2N &2N + ^ \®tV 

3=0 

Then a linear combination of the estimates (I4.99|) of Proposition 14.12l and the estimates (14.601) 
(14.611) of Proposition 14.91 gives 

rt rt 

V 2N ds <C £ £% N (0) + J~ 2 n{ 0) + C £ I \1&2N (^ 2 N T J~ 2 n) ds 

Jo 

+ £ [ {£2 N + JJ~ 2 n) ds + C £ I ll^-II^TV-i /2 ds (4.119) 

Jo Jo 

for any e > 0 and a corresponding constant C £ > 0. 

Next, we show that £% N is comparable to £f N and that comparable to T> 2 N■ We begin 

with the result for the energy. 

Proposition 4.13. Let £f N and £% N be as defined in (14.11) and (14.1181) respectively. It holds 
that 

£ 2 N < £ 2 N + £%. N (£ 2 N + ^2n) ■ (4.120) 


Proof. We compactly write 

X 2N = ||v 47V - 2 g 1 h; + ||v 4Ar_2 G 2 ||g + ||v 47V - 2 G 3 llL + llvf^G 4 " 2 


11/2 


1/2 


(4.121) 


We first estimate dfu for j = 0,... , 2N — 1. The key is to use the elliptic regularity theory 
of the following two-phase Lame system derived from (14.28|) : 


—p,Au — (fj/3 + //) V div u = G 2 — pdtu — pS7 ( h'(p)q ) in 
-§(-u)e 3 = ( -P'{p)q + pigrf+ - a + A if r] + )e 3 + G on S + 
- |S(u)] e 3 = (lP'(p)qj + [p] grf- + cm A*r/_)e 3 - G 3 _ on S_ 
[u| =0 on £_ 

U- = 0 on Sfe. 


(4.122) 


We let j = 0,..., 2 N — 1 and then apply <9^ to the problem (14.1221) and use the elliptic estimates 
of Lemma lA.51 with r = 4 N — 2 j > 2, by (14.1211) and the trace theory to obtain 

2 

//4JV-2j-3/2( S ) 


d{u 

2 

< 

dfG 2 

2 

+ 

d{ +1 u 

2 

+ 

d{g 

2 

+ 

dig 


4N-2j ~ 


41V—2/—2 


41V—2J-2 


41V—2/ —1 



+ 


&tV 


< 


di +l i 


4IV-2/-3/2 
2 

41V-2(j+l) 


+ <T* 


dlv 


4IV—2/+1/2 


+ 


&>o 3 


+ £ 2 IV + %2N ■ 


4A-2/-3/2 


Using a simple induction based on the estimate (14.1231) . utilizing the ||<9 2JV 'u||q estimate contained 
in £f N for the base case, we easily deduce that for j = 0,..., 21V, 

~ £2N + %2N- 


(4.123) 


diu 


41V—2/ 


(4.124) 


We then estimate Ofy and dP t r\ for j = 1 ,..., 2 N to get an improvement. By the first equation 
of (14.281) . using the estimates (14.1241) and ([4.1211) . we have that for j = 1,... , 2 N, 


dig 


< 


4N-2j+l 


di S 


df X u 


41V-2/+2 
2 


+ 


41V— 2{j- 1) 


+ 


di^G 1 


di^G 1 


41V-2/+1 

2 

41V—2(j—1)—1 


£ 2 N + ■ 


(4.125) 
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Now by the kinematic boundary condition 

d f r) = u 3 + G 4 on £, 

we have that for j = 1,..., 2 N, by the trace theory, (14. 124|i and (|4.121l) . 


(4.126) 


di: r l 

2 

< 

^ 1 u 3 


47V-2/+3/2 



< 

£>i~ i u 

4 


#4iV-2j+3/2( E ) 


+ 


+ 


&l~ l G 4 


-1^4 
2 


% G 


4AT-2/+3/2 


< 


&2N X‘2n■ (4.127) 


4JV-2(j-l)-l/2 

Summing the estimates (14.12411 . (14.12511 and (14.12711 . we conclude that 

&2N iS 4V + ■ (4.128) 

Using the estimate (14.3711 of Lemma [4.41 to bound (^ 2 IV + P2n), we then obtain 

flUHD from (144281) . " □ 

Next we consider a similar result for the dissipation. 

Proposition 4.14. Let T> 2 N an( ^ ^ 2 at fre as defined in (14.211 and (14.9811 respectively. It holds 
that 

T^2N 1$ T ) 2N + &2N (732N + &2N + -T^aO . (4.129) 

Proof. We compactly write 


y 2 iv = || + ||V 47V - 2 9 t G i || n + IIV^^G' 


+ llvf iV - 1 G 4 ||5 /2 + II v 4JV - 2 ^g 4 h; / 2 + u 2 ||V 4JV G‘ 


0 1 II v * "111/2 

- 4J V-2«-4||2 /2 + a 2|| v 4Ar G 4||2 /2 


(4.130) 


We now estimate the remaining parts of T> 2 N not contained in T>f N . We divide the proof into 
several steps. 

Step 1 - djq estimates 

We first notice that by the first equation of (I4.28p . 


— IMlLv + ll^* 1 11 4A/"_1 ~ ^ 2N + ^2N, 


(4.131) 


and for 2 < j < 2N + 1, 

d{ q 


VI 

CM 

d ] t 1 u 

2 

+ 

di^G 1 

4N-2j+2 


4AT-2/+3 


< 

di~ 4 u 

2 

+ 

di^G 1 



4AT-2/+3 



4N—2j+2 
2 

4N-2j+2 


< T>2N + y2N- 


(4.132) 


Step 2 - dfr) estimates 

We now derive estimates for time derivatives of q. For the term d° t p for j > 2 we use the 
kinematic boundary condition 

dtjj = u 3 + G 4 on £. (4.133) 

Indeed, for j = 2,..., 2 N + 1 trace theory and (14.1301) imply that 

2 


dip 

2 

< 

di 1 U 3 

2 

+ 

di~ l G A 


4N-2j+b/2 ~ 


H 4JV-2j+5/2( S ) 



4N—2j+5/2 


< 


di 1 u 

2 

+ 

^'- 1 G 4 


4AI—2(j—1)+1 



< 


P*2N + 3^2 A ■ (4.134) 


4AI—2(j—l)+l/2 

For the term dtp, we again use (I4.133I1 . trace theory, and (|4.130l) to find 

ll^ ? ?ll4Af+l/2 + ll^i r /ll4Ar-l/2 ~ + °' 2 ) ll w 3 ||^4JV+i/2( S ) + U 2 ||G 4 || 47V+1 /2 + 11 114A^—1/2 


~ IMLat+1 + y2N ^ £*2 N + y2N- 

Step 3 - V*?? estimates 

In this step we use instead the dynamic boundary condition 

- cr + A*r/+ + piOT+ = P' + {pi)q + - 2 p + d 3 u 3}+ - p' + di vu + - G : { + on £ + 


(4.135) 


(4.136) 
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and 

- cm A*p_ = - [p'(p)g] + 2 I/i( 9 3 « 3 l + [p' div w] - G| _ + [p] on X_. (4.137) 

Notice that at this point we do not have any bound of q on the boundary X, but we have 
bounded V(h'(p)q) in Ll. As such, we first apply V* to (14.1361) and (14.1371) and then we employ 
the standard elliptic theory. This, trace theory, and (14.1301) then provide the estimate 


U 2 || V *? 7|| 4W - +1 / 2 + 11 11 4 JV— 


3/2 


II /II" 2 || o n2 2 

||v*(/l (p)q) ||^ 4 iV- 3 / 2 ( S ) + ||V*Vu||^ 4 JV- 3 / 2 ( S ) + 11 ^3 11 4 JV- 3/2 + *V- II 4 JV— 3/2 

~ ||V(/l , (p)<?)||4 W -_ 1 + |Ml4N- +1 + ||G 3 


I 4 IV—1/2 II 7 ?- II 4 JV—1/2 


< X>2A + 3^2iv + HP-H 4 JV- 1/2 • (4.138) 

Consequently, summing the estimates (14.1311) . (I4.132|) . (|4.134|) . (14.1351) and (I4.138|) . we con¬ 
clude 

^2 N ~ ^2 N + y 2 N + \\v~ 11 47 V— 1 /2 • (4.139) 

Using the estimate (|4.38l) of Lemma f4.4l to bound T 2 V 1$ £%n (^ 2 N + J~ 2 n), we then obtain 
dUMD from (144391) . " □ 

Finally we are ready to prove Theorem 14.11 

Proof of Theorem \f.l\ By the estimates (14.1201) of Proposition 14.131 and the estimates (14.1291) 
of Proposition 14.141 we can improve the inequality (14.1191) to be 

&2Nif) + J~'2N(t) + f Df^rds TC £ £f N (0) + P2n(^) + C £ f yj £% n (fDfN + ^2N + J~2n) ds 

Jo Jo 


+ s 


[ (£ 2 n JF 2 n) ds + C £ 
Jo 


I4JV—1/2 


ds. 


Sobolev interpolation on X allows us to bound 


C E 


4JV-1/2 — 


< e 77 J 


4N 


+ C £ ||r /_||5 < e£fjv + C £ ||r /_" 2 


0 • 


We can thus refine the inequality ()4. 1 191) to be (14.51) . 

5. Nonlinear instability 


(4.140) 

(4.141) 

□ 


5.1. Restated estimates. In the following, we take A = A defined by ([2.341) when cr_ > 0, 
while we take A = A* defined by (12.351) when cm = 0. In each case, we have that 4 < A < A. 


We define the norm 


II 00 


appearing in Theorem 11.21 by 


\l(Q,u,v)\ioo : = V^ 2 N + ^n (5.1) 

for an integer N > 3, where Sf N and JF 2 N are given by (14.11) and (14.31) . For notational conve¬ 
nience, we denote 

U := (q,u, rf). (5.2) 

We now restate the main results of the previous Sections in our new notation. 

Proposition 5.1. Let the norm I-Iqo be given by (EU). Then we have the following. 

(1) There is a growing mode U* := (q*,u*,rj*) satisfying ||r/* || f) = 1, |||f7*||| 00 = Ci < 00 , 
and e xt U* is the solution to m- 

(2) Suppose that U(t) is the solution to (I4.28|) . There exists a small constant S such that if 
I Uiff) Iqq < 5 for all t € [0, T], then there exists C$ > 0 so that the following inequality 
holds for t € [0, T]: 

l|t/(t)lllm < Cj ||t/(0)|£, + 1 [ It/WIJ, is 


+ Cs 


fm 

Jo 


II00 


ds T C $ 


[ ll ? 7-(s)llo ds - 
Jo 


(5.3) 
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(3) There exists C 2 > 0 so that 

I \V-(t) - te x W-\\o < C 2 e M |||C/(0) - tti*||| 00 + C 2 f ||t/(s)l»oo ds 

J 0 

+ C *Jf* e2A(i_s) I^WlSo 1^) - ^ XtU %o ds • M 

Proof. Statement 1 follows from Theorem 12.71 Statement 2 follows from the estimates (|4.5p 
of Theorem 14.11 by taking e = A/2 and then taking 5 sufficiently small to absorb the term 
yj£f N V 2N 011 the r ight hand side by the dissipation. 

We now prove Statement 3 by using Theorem 13.31 We observe that U(t) — ie xt U* solve the 
problem (13.11) with initial data U( 0) — lU * and the force terms G l given by ()4.29f) (14.36[) . Then 
Statement 3 follows from (13.201) by noticing that 

= ||G' 1 ||i + \\d t G 2 \\ 2 0 + H^^Ho + l|G' 4 || 2 < iUfoo ■ ( 5 - 5 ) 

We thus conclude the proposition. □ 


5.2. Local well-posedness. Thus far we have not elaborated on the local well-posedness the¬ 
ory for our problem that we developed in our companion paper [10]. In the result we will refer 
to the “necessary compatibility conditions” required for the local well-posedness in our energy 
spaces. These are cumbersome to write out explicitly, and we refer to [TO] for the explicit 
statement. 

Theorem 5.2. Suppose that the initial data U( 0) satisfies the necessary compatibility condi¬ 
tions. There exist 5 q,T > 0 so that if 

r(o)ioo < £0, 

then there exists a unique solution U(t ) to (I4.28P on [0, T] that satisfies the estimate 

l^)loo^vT+r|^(0)IL 

for all t E [0, T], 

Proof. The theorem can be deduced readily from Theorem 2.1 of nnj. Indeed, Theorem 2.1 
of PE] is stated in more general form, where we only require ||r/o|| 4 jv-i /2 to be small and no 
smallness condition is imposed on uq or q 0 . We record this version of local well-posedness so 
that it can be adapted directly in our instability analysis. □ 


(5.6) 

(5.7) 


5.3. Data analysis. In order to prove our nonlinear instability result, we want to use the linear 
growing mode in Proposition [37T] to construct small initial data for the nonlinear problem (I4.28p . 
Since we are involved in the higher-order regularity context, we cannot simply set the initial 
data for the nonlinear problem to be a small constant times the linear growing modes. The 
reason for this is that the initial data for the nonlinear problem must satisfy certain nonlinear 
compatibility conditions in order for us to guarantee local existence in the space corresponding 
to norm |||-lloo’ which the linear growing mode solutions do not satisfy. 

To get around this obstacle, we note that the nonlinear problem is slightly perturbed from 
the linearized problem and so their compatibility conditions for the small initial data should 
be close to each other. We are able to produce a curve of small initial data satisfying the 
compatibility conditions for the nonlinear problem which are close to the linear growing modes. 


Proposition 5.3. Let U* be the linear growing mode stated in Proposition I5.il Then there 
exists a number to > 0 and a family of initial data 

U t l = lU* + t 2 t/(t) (5.8) 

for l E [0, to) so that the followings hold. 

1. Uq satisfy the nonlinear compatibility conditions required by Theorem, 1 5 . 21 for a solution to 
the nonlinear problem (14.281) to exist in the norm ||| - |q 0 - 
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2. There exist C 3 ,C 4 > 0 independent of i so that 


U(i) < c 3 

00 

(5.9) 

and 


Moo < Cal 2 - 

(5.10) 

Proof. See the abstract argument before Lemma 5.3 of [9]. 

□ 


5.4. Proof of Theorem ll.2t With Propositions 15. II Theorem [572] and Proposition [5J3] in hand, 
we can now present the 


Proof of Theorem \1.2[ Recall the notation (15.21) . First, we restrict to 0 < l < lq < 9q, where lq 
is as small as in Proposition 15.31 and the value of 6q is sufficiently small to be determined later. 
For 0 < l < to, we let Uq be the initial data given in Proposition 15.31 By further restricting l 
we may use (15.101) to verify that (15.61) holds, which then allows us to use Theorem 15.21 to find 
U L (t), solutions to the system (|4.28l) with 

U L \ t=0 = U& = CU * + (5.11) 

Fix 5 > 0 as small as in Proposition 15.11 and let Cs > 0 be the constant appearing in 
Proposition 15.II for this fixed choice of <5. We then define 5 = min{<5, Denote 


T* = sup 


{* = 


117 ' 


oo < <5, for 0 < t < s 


and 







T 

** = sup js : r?7(f) 

0 < 2 te x \ 

for 0 

VI 

With to small enough, (|5, 

.101) and (|5.7|) guarantee that T 

* and 

p>** 

bv (fL35l). 

Then for all t 

< min {T L ,T*,T**}, ■ 

we deduce 

! from 

the 

15.11 the definitions of T* 

and T**, and (|5TTJD 

that 




iircoioo 

< C S |(7oloo + ^j‘ 

0 to 

0 

ds 




+ 

0 00 

0 

ds + Cs J 

r V- 

0 

(s) 


(5.12) 

(5.13) 


< 


< A 


+ SC s 


f lU\a) 

Jo 


f 

Jo 


ll2 J_ , ^ ri .2 , 2\t 


(5.14) 


U\s)\l z 00 ds + C S C 4 C + 


2A 


loo ds + C 5 i 2 e 2Xt . 


for some constant C 5 > 0 independent of 1 . We may view (15.141) as a differential inequality. 
Then Gronwall’s lemma implies that 

rt 


III^(t) 1 00 < C 5 i 2 e 2Xt + C 5 i 2 e xt [ \e Xs ds 

Jo 

< C b t 2 e 2Xt + C b i 2 e 2Xt = 2 C 5 c 2 e 2 
We then deduce from Proposition 15.11 and (15.151) that 

rt 


(5.15) 


_(t) - ie xt rf_ ||o < C 2 e At l 2 U(l) 


00 


+ c 2 


ds 


l fL 


+ CW / e^t-s) |||J7 t (s)|||oo |C/ l ( 5 ) - ^U% 0 ds 


< C 2 C 3 e A h 2 + C 2 I 2 C 5 Pe z * s + C 2 


/ 

7o 


,2„2As 


e 2A J-s)2C 5 t 2 e 2Xs (^2C5te Xs + Cne Xs )ds 
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Here we have used the fact that A < 2A. 

Now we claim that 

T L = min{T\ T*, T**} (5.17) 

by fixing 8 q small enough, namely, setting 

e ° =min {vWsri 4 }' (5 ' 18) 

Indeed, if T* = min {T L ,T* ,T**}, then by (15.151) . we have 


l^cniloo < y/2C 5 ^ < V^C~5ie XTL = y/2Cle 0 < 6 - < 6, (5-19) 

which contradicts to the definition of T*. If T** = min{T t , T*, T**}, then by (15.161) and the 
fact that A/2 < A < A, we have that 

\\V-(T **)\\ 0 < te AT " b -\\ 0 + I \V-(T**) ~ ie XT **rf _|| 0 

< ie AT " \U_\\+2C 6 L 2 e 2XT ** + Ceiled 

0 , (5.20) 

< ie XT (1 + 2<7 6 te AT ‘ + C^ XTL ) 

< te XT ** (1 + 2 C 6 8 0 + Cgv^) < 2 ie AT “, 

which contradicts to the definition of T**. Hence (15.171) must hold, proving the claim. 

Now we use (15.161) again to find that 

||rtcn || 0 > ,e AT ‘ b *_|| 0 - I \r,L(r) - || 0 

> te AT ‘ — 2 C G L 2 e 2XTL - C§L^e^ xt 

>e 0 - 2 c 6 e 2 - c 6 e$ > y. 

This completes the proof of Theorem 11.21 


(5.21) 

□ 


Appendix A. Analytic tools 


A.l. Poisson extensions. We will now define the appropriate Poisson integrals that allow us 
to extend rj±, defined on the surfaces to functions defined on H, with “good” boundedness. 

Suppose that = T 2 x {j}. where T 2 := ( 2ttL\T ) x ( 2 ttL 2 T). We define the Poisson integral 
in T 2 x ( oo, j) by 

5e(A Z)x(^ Z) 


where for £ G (L 1 1 Z) x (L 2 1 Z) we have written 



fix') 


e -i£-x' 
27T y/ L 1 L 2 


dx'. 


(A.2) 


Here ” stands for extending downward and “j” stands for extending at X 3 = j, etc. It is 
well-known that V- t j : H S (Y, + ) —> H s+1 / 2 ( T 2 x (—00, j)) is a bounded linear operator for s > 0. 
However, if restricted to the domain H, we can have the following improvements. 


Lemma A.l. Let V— jf be the Poisson integral of a function f that is either in H q (Y>+) or 
iT , ~ 1 / 2 (E + ) for q G N = {0,1,2,...}, where we have written H S (E + ) for the homogeneous 
Sobolev space of order s. Then 

HV^-j/llo < \\f\\% q - 1/2{T2) and \\V q V-jf \\ 0 < \\f \\% q{T2) . (A.3) 

Proof. See Lemma A.3 of [7j. □ 
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We extend rj+ to be defined on Ll by 

fj + (x',x 3 ) = V+rj + (x',x 3 ) := V- }j r] + (x', x 3 ), for x 3 < j. 


(A.4) 


Then Lemma I A. II implies in particular that if t] + G H s 1 / 2 (S + ) for s > 0, then r] + G H s (Ll). 
Similarly, for E„ = T 2 x {0} we define the Poisson integral in T 2 x (—oo,0) by 


V-,of(x) = 


E 


J£-x' 




‘1'K\J L\Li 


m- 


(A.5) 


It is clear that V—o has the same regularity properties as V—j. This allows us to extend to 
be defined on fl_. However, we do not extend ??_ to the upper domain by the reflection 
since this will result in the discontinuity of the partial derivatives in x 3 of the extension. For 
our purposes, we instead to do the extension through the following. Let 0 < Ao < Ai < • • • < 
A m < oo for m G N and define the (m + 1) x (m + 1) Vandermonde matrix V(Ao, Ai,..., A m ) 
by V(Ao, Ai,..., \ m )ij = (—A j ) 1 for i,j = 0 , It is well-known that the Vandermonde 

matrices are invertible, so we are free to let a = (ao, or, •.., a m ) T be the solution to 

V(Ao, Ai,..., A m ) a — Qmi Qm — (1) 1) • • • j 1) ■ (A. 6 ) 

Now we define the specialized Poisson integral in T 2 x (0, oo) by 


'P+,of(x) 


E 

teiL^zML^z) 


e i^x' 

27 T\/ L1L2 


m 


3=0 


(A.7) 


It is easy to check that, due to (IA. 6 I) . d l 3 V+fif(x',0) = d l 3 'P^.of(x',0) for all 0 < l < m and 
hence 

d a V+fif{x', 0) = d a V- t of(x', 0), Va G N 3 with 0 < |a| < m. (A. 8 ) 

These facts allow us to extend ?]_ to be defined on H by 


rj-(x r , x 3 ) = V-rj-{x',x 3 ) 


'P+, oV-(x',x 3 ), x 3 >0 

V- fi V-(x',x 3 ), x 3 <0. 


(A.9) 


It is clear now that if f?_ € H s ~ l / 2 ( E_) for 0 < s < m, then r/_ G H S (Q). Since we will only 
work with s lying in a finite interval, we may assume that m is sufficiently large in (IA. 6 I) for 
fj- € H s (fl) for all s in the interval. 


A.2. Estimates of Sobolev norms. We will need some estimates of the product of functions 
in Sobolev spaces. 

Lemma A.2. Let U denote a domain either of the form or of the form E±. 

(1) Let 0 < r < si < S 2 be such that s\ > n/2. Let f € H Sl (U), g € H S 2 (U). Then 
fg E H r (U) and 

\\fg\\ H r ~ ll/ll/pu ■ (a. 10 ) 

(2) Let 0 < r < s\ < S 2 be such that S 2 > r + n/2. Let f G H S 1 (U), g G H S 2 (U). Then 
fg G H r (U ) and 

ll/5 , lln r ~ ll/ll// s i lbll// s 2 ■ (A-ll) 

Proof. These results are standard and may be derived, for example, by use of the Fourier 
characterization of the H s spaces and extensions. □ 
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A.3. Coefficient estimates. Here we are concerned with how the size of if can control the 
“geometric” terms that appear in the equations. 

Lemma A.3. There exists a universal 0 < 6 < 1 so that if ||^7]|§/2 — ^ then 

1 


\J ~ 1 llL°°(n) + 


L°°(G) + 1111 L°° (P) — 2 ’ 

1 


r) A ||A l|li°o(r) — ^ ant ^ 
< 1 . 


II-AA- 1| 

\\ k \\ l ™( q .) f ii^iiL°°(n) 

Also, the map 0 defined by (11.181) is a diffeomorphism. 

Proof. The estimate (1A.12D is guaranteed by Lemma 2.4 of [7]. 


(A.12) 


□ 


A.4. Korn inequality. Consider the following operators acting on functions u: 


Du = Vu + Vu r and D°u = Du — 


2 divu 

3 


I, 


or in components 


D Uij = diUj + djUi and 


u. 


v 


— diUj T djUi 




Jiy 


Note that tr(D°u) = 0. As such, the operator D° is referred to as the “deviatoric part of the 
symmetric gradient.” 

Now we record a version of Korn’s inequality involving only the deviatoric part, D°, that we 
will use for layered domains f2±. 


Proposition A.4. There exists a constant C > 0 so that 

IH|?<C||D°u||o (A.13) 

for all u E with [u] = 0 along £ and u_ = 0 on £&. 

Proof. We refer to Proposition A.8 of m- □ 


A.5. Elliptic estimates. Here we consider the two-phase elliptic problem 


—fjAu — (p,/ 3 + //)V div u = F 2 

in D 

-§(u + )e 3 = F'l 

on £ + 

- [S(u)] e 3 = -Fl 

on £_ 

M = o 

on £_ 

u_ = 0 

on 


(A.14) 


We have the following elliptic regularity result. 

Lemma A. 5. Let r > 2. If F 2 € H r ~ 2 (Ll),F 3 € Lf r-3 / 2 (£) ; then the problem (1 A. 1411 admits a 
unique strong solution u E H r (Ll). Moreover, 

IMIr ~ ll i7l2 |lr-2 + 11^11,-3/2 ' (A.15) 

Proof. We refer to m Theorem 3.1] for the case of two-phase Stokes problem, but the proof is 
the same here. It follows by making use of the flatness of the boundaries S-t and applying the 
standard classical one-phase elliptic theory with Dirichlet boundary condition. □ 


We let G denote a horizontal periodic slab with its boundary dG (not necessarily flat) consist¬ 
ing of two smooth pieces. We shall recall the classical regularity theory for the Stokes problem 
with Dirichlet boundary conditions on dG, 

—frAu + Vp = / in G 

< di vu = h in G (A. 16) 

u = tp on dG. 

The following records the regularity theory for this problem. 
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Lemma A.6. Let r > 2. If f € H r 2 (G),h € H r 1 (G),ip € H r 1//2 (i dG) be given such that 


then there exists unique u € H r (G),p 
IMI H r (G) + l|Vp||j ? r-2( G ) < | 
Proof. See El. 


/* = / ■ 

iG J dG 


P ■ V, 


(A.17) 


€ H r 1 (G)(up to constants) solving (1A.16I) . Moreover, 

/II JR— 2 (G) + IWliP-RG) + IMIiR—V^dG) • (A.18) 

□ 
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